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que podem interagir com a 2PAM e a distâncias em A. (B) Distâncias 
destacadas em (A), calculadas ao longo da simulação 99 
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HuAChE-GA-DZP(C), onde estão destacados a DZP e os resíduos mais 
próximos dela na saída do sítio ativo com o SGA no fundo do poço 100 

FIG. 3.33 Ultimo quadro referente ao complexo HuAChE-GA-DZPanc(C) com o 
sítio ativo representado em superfície de van der Waals , a DZPanc e o 
SGA em varetas. Onde em (A) tem-se a vista lateral e em (B), a vista 


superior 100 

FIG. 3.34 Cálculo da variação de energia potencial eletrostática para (A) 2PAM, 
(B) DZP, (C) DZPanc e (D) sítio ativo da HuAChE 101 
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RESUMO 


O projeto de reativadores da acetilcolinesterase é fundamental para o 
desenvolvimento de antídotos contra os agentes neurotóxicos de guerra química, 
assim como o uso dessa enzima em equipamentos de defesa química. Entretanto, o 
mecanismo de reações de reativação e a característica estrutural que determina a 
eficiência dos reativadores são pouco conhecidos. Neste trabalho, fizemos estudos 
por dinâmica molecular do cristal da acetilcolinesterase humana (HuAChE) inibida 
pelo tabun, complexado com o antídoto pralidoxima (2PAM) e com seus deazo 
análogo nas formas neutra e aniônica. Os resultados mostraram que a carga 
positiva da 2PAM é crucial para sua permanência e para sua entrada no sítio ativo 
da HuAChE enquanto a deazapralidoxima (DZP) neutra e aniônica são repelidas 
quando próximas à entrada. Os resultados e análises por dinâmica da 2PAM, dentro 
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e fora do sítio ativo da enzima, podem ser úteis para o desenvolvimento de novos 
antídotos potentes contra os agentes neurotóxicos. 


ABSTRACT 


The design of efficient acetylcholinesterase (AChE) reactivators is fundamental 
for the development of antidotes against neurotoxic Chemical warfare agents, as well 
as for the use of that enzyme in Chemical defense equipment. However, the 
mechanism of the reactivation reaction and the structural characteristics that 
determine the efficiency of a reactivators are poorly understood. In this work, we 
carried out a molecular dynamics (MD) study of the crystallographic human 
acetylcholinesterase (HuAChE) inhibited by tabun in complex with the known 
antidote pralidoxime (2PAM) and with its deaza analogue in the neutral and anionic 
forms. Results show that the positive charge of 2PAM is crucial for its permanence 
and its entrance inside HuAChE active site while the uncharged and anionic 
deazapralidoxime (DZP) moves away from the entrance. Those results and the 
analysis of the dynamic behavior of 2PAM inside and outside the enzyme could be 
very useful for the design of new potential antidotes to neurotoxic agents. 
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1 INTRODUÇÃO 

1.1 POSICIONAMENTO E ABORDAGEM DO PROBLEMA 

Uma das deficiências atuais da defesa química reside na limitada eficácia dos 
antídotos disponíveis contra os agentes dos nervos, a classe de agentes químicos 
que possui maior toxidez. No desenvolvimento deste trabalho, aplicaram-se as 
técnicas de Modelagem e Dinâmica Molecular (DM) para um estudo das interações 
e da trajetória entre um dos antídotos mais usados, a pralidoxima, e seu análogo 
protonado e desprotonado com os resíduos do sítio ativo da enzima, onde o 
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caminho percorrido pelo fármaco e os similares têm como referência à posição do 
aminoácido fosforilado pelo agente Tabun. 

Este trabalho faz parte de um projeto maior que tem como objetivo final o 
desenvolvimento de novos e mais eficientes antídotos para a intoxicação com 
compostos organofosforados neurotóxicos. Este projeto envolve um grande número 
de etapas, algumas das quais já foram alcançadas. Estas etapas incluem, desde a 
obtenção de um modelo virtual tridimensional da acetilcolinesterase humana e da 
enzima inibida com os diferentes organofosforados, ao desenvolvimento de um 
modelo para simulação da reação de desfosforilação da enzima com oximas, 
levando-se a estudos da dinâmica e do atracamento molecular das oximas 
atualmente usadas como antídotos. Assim, pelo uso de todas as informações 
obtidas para o projeto racional de novos potenciais antídotos, busca-se a síntese e a 
avaliação biológica dos mais promissores. 

1.2 OS AGENTES DE GUERRA QUÍMICA 

Desde a antiguidade, formas primitivas do emprego de compostos químicos 
como agentes de guerra química, devido a seus efeitos tóxicos no homem, animais 
ou plantas, é relatado na literatura (WHO, 1970). Idéias sugeridas, mas não 
colocadas em prática, incluem, por exemplo, o uso de munições com cianeto pelos 
ingleses na Guerra do Criméia e munições com cloro na Guerra Civil Americana. 
Posteriormente, com o desenvolvimento da indústria química do século XX, a 
utilização destes agentes tóxicos em guerra química tornou-se mais ampla. Em 
função disso, surgiu um ramo da defesa militar e civil, a chamada DQBN (Defesa 
Química, Biológica e Nuclear), que tem se preocupado com o desenvolvimento de 
sistemas de defesa que minimizem os efeitos tóxicos dos agentes químicos de 
guerra, biológicos e nucleares, na capacidade operacional dos exércitos e na 
população exposta (RICKETT, 1987). 

O início da guerra química moderna foi em 1915, na Primeira Guerra Mundial, 
em que os alemães utilizaram gás cloro contra tropas aliadas em Ypress, Bélgica, 
levando assim ao surgimento das primeiras máscaras contra gases e, 
conseqüentemente, à necessidade de agentes mais tóxicos (como o fosgênio) ou 
cujas características físico-químicas dificultassem sua retenção pelas máscaras 
(como a cloropicrina). Outra evolução foi o desenvolvimento de agentes que causam 
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danos à pele (como as mostardas e os vesicantes arsenicais). Ao final da guerra, a 
utilização de agentes químicos por ambos os lados causou aproximadamente 100 
mil mortos e 1,3 milhões de feridos (FOA, 1992). 

Entre a década de 30 e a II Guerra Mundial, uma nova classe de agentes de 
guerra foi sintetizada na Alemanha. Estes, denominados neurotóxicos, por agirem 
sobre o sistema nervoso, tinham uma toxidez muito maior que os conhecidos 
anteriormente. Em 1950, novos agentes desta classe foram sintetizados na 
Inglaterra. Estes possuem uma grande facilidade de síntese. Em virtude disso, a 
ONU não tem logrado êxito no controle da produção destes agentes químicos de 
guerra, o que favorece também problemas na área de defesa Civil. 

A defesa contra um ataque por agentes químicos é realizada em três níveis. O 
primeiro, e melhor, é evitar o contato ou permanecer em uma área não exposta a 
contaminação. O segundo é o uso de equipamentos de proteção individual e 
coletiva para evitar o contato direto com os agentes e a conseqüente intoxicação. 
Em último caso, vários estágios de tratamento médico especializado (pré- 
tratamento, antídotos auto-administráveis etc.) buscam minimizar os efeitos da 
intoxicação. 

Atualmente, a maior preocupação em relação à guerra química está focalizada 
nos agentes neurotóxicos. Estes compostos agem principalmente na 
acetilcolinesterase (AChE), enzima responsável pelos processos de transmissão 
dos impulsos nervosos. 

1.2.1 Os Agentes Neurotóxicos 

1.2. 1.1 Histórico 

São os representantes mais tóxicos do grupo de substâncias anti-ChE inibidoras 
da enzima acetilcolinesterase, incluindo também inseticidas organofosforados, 
carbamatos, haletos de sulfonila, proteínas naturais do grupo fasciculina e outros 
inibidores sintéticos reversíveis altamente específicos. A pesquisa dos efeitos da 
exposição a substâncias anti-ChE é associada principalmente à defesa química 
contra agentes neurotóxicos e ao tratamento da exposição acidental a inseticidas 
organofosforados, estimada em um milhão de casos por ano, principalmente em 
países do “Terceiro Mundo” (WHO, 1990). 
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O primeiro agente neurotóxico sintetizado foi o tetraetilpirofosfato (TEPP), em 
1854, mas no geral, o desenvolvimento de organofosforados, como agentes de 
guerra, só ocorreu oito décadas depois (SIDELL, 1992). 

Na década de 1930, a necessidade de aumentar a produção de alimentos 
estimulou a pesquisa de novos compostos com ação inseticida. Em particular, o 
cientista alemão Gerhard Shrader do laboratório I.G. Farben sintetizou em 1936 o 
N,N-dimetilfosforamidocianidato de o-etila (conhecido como Tabun), que 
apresentava uma toxidez aos mamíferos alta demais para permitir o seu emprego 
comercial como pesticida. O potencial militar do Tabun foi reconhecido pelos 
militares alemães que passaram a produzi-lo em grandes quantidades para o 
emprego em munições químicas, com a produção chegando a um total estimado de 
10.000 - 20.000 toneladas/ano ao final da guerra. Dois anos depois, Shrader 
sintetizou o metil fosfonofluoridato de isopropila (conhecido como Sarin), ainda mais 
tóxico que o Tabun. O interesse pela síntese destes agentes foi tão grande que em 
1944 um outro grupo de pesquisadores sintetizou o etil fosfano fluoridato de 
pinacolina (Soman). 

Na década de 50, mais um composto organofosforado foi desenvolvido pelos 
ingleses para uso como pesticida. Este composto, conhecido como Amilton, revelou- 
se tóxico demais para emprego comercial. A modificação do Amilton levou a 
descoberta do primeiro agente da chamada série V, o VX [S-(2-diisopropilaminoetil) 
metil fosfonotiodato de etila], Estes dois agentes apresentavam, ao contrário da 
série G (Tabun, Sarin e Soman), baixa volatilidade e uma toxidez por absorção 
cutânea quase igual à toxidez por inalação. Sarin e VX eram os agentes 
neurotóxicos padrões dos EUA, enquanto a União Soviética estocava Soman. O 
Tabun foi abandonado pelas principais potências em função de sua menor toxidez, 
mas sua maior facilidade de produção tornou-o atraente para países que pretendiam 
iniciar um programa de desenvolvimento de armas químicas. 

Hoje em dia, GA, GB, GD e VX (FIG.1.1.), são os quatro principais agentes 
neurotóxicos padronizados para uso em munições químicas e estocados em 
grandes quantidades por vários países (FRANZ, 2004). 
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V A 

FIG. 1.1 - Principais agentes neurotóxicos. 


No geral, estes agentes apresentam um átomo de fósforo quiral ligado a um 
radical RI (alquila ou dialquilamino), um radical R2 [alquila, cicloalquila, H(somente 
na série V), (CH 2 ) n N + R 3 (quando X=F)], formando dupla ligação com um grupamento 
Y (oxigênio ou raramente enxofre) e ligado a um grupo de saída X (como o flúor e o 
cianeto, por exemplo). 

O primeiro uso, comprovado destes agentes foi por volta de 1980, na guerra Irã- 
Iraque, onde as forças iraquianas empregaram o GA. Durante a Guerra do Golfo, 
tropas americanas foram expostas ao GB após destruírem um depósito iraquiano de 
munições químicas. Em 1994 e 1995, GB foi usado em ataques terroristas no 
Japão, dos quais o mais grave, no metrô de Tóquio, deixou 12 mortos e 5500 
feridos (TU, 1998; MASKAWA, 1995). 


1.2. 1.2 Principais Propriedades Físicas 


Apesar de serem conhecidos como “gases dos nervos” e de serem observados 
em forma de vapor, todos estes quatro agentes (GA, GB, GD e VX), apresentam-se 
sob a forma líquida em temperatura ambiente. O que os torna “gasosos” é 
meramente a arrebentação das munições que os contém ou o espargimento 
(SIDELL, 1997; BLACK, 1996). 

O VX é um agente de menor volatilidade (maior ponto de ebulição) o que o torna 
persistente no meio ambiente, podendo levar semanas para ser disperso. Já a série 
G (GA, GB e GD) é moderadamente volátil, dispersando no ambiente em poucas 
horas. Assim, apesar destes quatro agentes serem bastante tóxicos por todos os 



meios de absorção, a série G apresenta um maior risco por inalação, enquanto o VX 
apresenta maior risco por absorção cutânea. As propriedades físicas destes agentes 
são resumidas na TAB. 1.1. 


TAB. 1.1 - Propriedades físicas dos agentes neurotóxicos (SIDELL, 1997; BLACK, 1996). 


Propriedade 

GA 

GB 

GD 

VX 

Ponto de Ebulição(°C) 

230 

158 

198 

298 

Pressão de Vapor (mm Hg) (20°C) 

0,035 

2,1 

0,34 

0,0004 

Densidade 

Vapor (ar=1) 

5,6 

4,86 

6,3 

9,2 

Líquido (g/ml) (20°C) 

1,077 

1,10 

1,01 

1,013 

Viscosidade (cp) (20°C) 

2,77 

1,54 

- 

12,2 

Solubilidade em água (%) (20°C) 

9,8 

Miscível 

2,1 

3,0 


1.2. 1.3 Propriedades Químicas 


Os estudos acerca das propriedades químicas dos agentes dos nervos têm-se 
concentrado nas reações referentes à sua detecção, descontaminação, 
estabilização e comportamento em sistemas biológicos (BLACK, 1996). No caso do 
sarin e do soman, resume-se meramente a quebra da ligação P-F por ataques 
nucleofílicos. Para o tabun e o VX, a ligação a ser quebrada inicialmente varia com 
os reagentes e condições de reação. Para todos estes quatro agentes o nucleófilo é 
um resíduo serina presente no interior do poço do sítio ativo da acetilcolinesterase, 
sendo este a base de sua toxidez (FRANZ WOREK, 2004). 

Valores letais para a contaminação a partir de organofosforados para várias 
espécies de seres vivos encontram-se disponíveis na literatura (MARRS, 1996). 
Porém, para os seres humanos estes valores são estimados. Essas estimativas são 
apresentadas na TAB. 1.2, onde o ácido cianídrico foi tomado como referência e a 
LD 5 o é a dose capaz de matar 50% da população exposta. No caso de exposição ao 
vapor, a dose é dada pelo produto C x t (concentração do agente no ar x tempo de 
exposição) (SIDELL, 1997; THE WEDNESDAY REPORT, 2004). 


TAB. 1.2 - Toxidez estimada dos agentes neurotóxicos ao homem 


Agente 

LCt 5 o (inalação) 

(mg. min/m 3 ) 

LDso (pele) 

(mg) 

Tabun (GA) 

400 

1000 

Sarin (GB) 

100 

1700 

Soman (GD) 

70 

50 

VX 

50 

10 

HCN 

2500-5000 

- 


*LD 5 q= dose capaz de matar 50% da população exposta, por absorção cutânea. 
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l_Ct 5 o= dose capaz de matar 50% da população exposta, por via respiratória. 


1.3 AACETILCOLINESTERASE 

As colinesterases definem uma família de serina hidrolases, dentro da 
superfamília de alfa/beta (a/p) hidrolases (OLLIS, 1992; NARDINI, 1999). As 
colinesterases são divididas em acetilcolinesterases (AChE) e butirilcolinesterases 
(BuChE), que apresentam diferentes especificidades para seus substratos naturais 
(ésteres de colina) e inibidores (TAYLOR, 1994). 

A acetilcolinesterase (AChE) é uma enzima responsável pela terminação da 
transmissão dos impulsos nervosos nas sinapses colinérgicas através de uma 
rápida hidrólise do neurotransmissor acetilcolina (ACh). E é por esta sua função 
chave, que ela se torna um dos alvos a serem usados em contaminações por gases 
neurotóxicos, pesticidas, venenos de cobras e também no tratamento com fármacos 
anti- Parkinsons. 

Existe uma grande variedade de neurotransmissores no sistema nervoso central 
(SNC). Dentre eles, os dois mais importantes, e que participam dos processos do 
sistema nervoso periférico, são a acetilcolina e a noradrenalina, ver FIG.1 .2. 


O 

C ^ Wlci 
H 3 c o 

Acetilcolina 

FIG.1.2 - Os dois principais neurotransmissores do sistema nervoso periférico. 

Ações referentes ao sistema nervoso periférico podem ser divididas em dois 
sistemas: somático e autônomo. O sistema autônomo pode ser classificado como 
simpático ou parassimpático. No sistema somático as estimulações se concentram 
nos músculos esqueléticos. Já no autônomo, tem-se a noradrenalina, pertencente à 
classificação autônoma simpática, que se concentra nos órgãos e principalmente 
nos músculos cardíacos, e a acetilcolina, pertencente ao sistema autônomo 
parassimpático, que, diferente da noradrenalina, reage com receptores específicos 
(PATRICK, 2002). 


HO 



R = Me Adrenalina 
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1.3.1 O Sistema Nervoso Periférico 


Como o nome indica, é uma parte do sistema nervoso que está na periferia do 
sistema nervoso central (SNC, cérebro e coluna vertebral). 

Existem algumas divisões e subdivisões do sistema periférico. As primeiras 
distinções que se podem fazer são as seguintes: 

• Nervos sensoriais (nervos que recebem mensagens no sentido do corpo para 
o SNC). 

• Nervos motores (nervos que levam a mensagem do SNC ao resto do corpo), 
que subdividem-se em: 

a) Sistema nervoso motor somático. 

b) Sistema nervoso motor autônomo. 

c) Sistema nervoso entérico. 

1.3. 1.1 Sistema Nervoso Motor Somático 

Não há sinapse e o neurotransmissor nas junções neuromusculares é a 
acetilcolina. Neste sistema, as mensagens são levadas do SNC para os músculos 
esqueléticos, onde o resultado final é a contração destes. 

1.3. 1.2 Sistema Nervoso Motor Autônomo 

Neste caso os nervos conduzem os impulsos no sentido do SNC para os 
músculos lisos, cardíacos e medula. Este sistema pode ser dividido em dois 
subgrupos (já comentado no primeiro parágrafo após a FIG.1.2): parassimpático e 
simpático. 

1.3. 1.3 Sistema Entérico 

Este terceiro subgrupo do sistema nervoso periférico está localizado nas paredes 
do intestino. Recebe mensagens tanto dos nervos autônomos simpáticos quanto 
dos parassimpáticos, além de responder a efeitos locais como reflexos importantes 
nas funções gastro intestinais (Gl). 
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Na FIG.1.3, é representado um esquema bastante sintético e explicativo dos três 
subsistemas do sistema nervoso periférico, sendo que os receptores nicotínicos são 
aqueles estimulados, por exemplo, pela nicotina (presente no tabaco) e os 
muscarínicos pela muscarina, encontrada nos cogumelos venenosos. 


SNC 

(Somático) 


ACh 

(N) 


Músculos 

esqueléticos 

SOMÁTICO 



(N) 



AUTONOMO 


Músculo Liso 
Músculo Cardíaco 


FIG.1 .3 - Nervos motores do sistema nervoso periférico. N, receptores nicotínicos; 
M, receptores muscarínicos (PATRICK, 2002). 


1.3.2 A Acetilcolina 


Como já foi comentada na seção 1.3, a acetilcolina tem grande importância na 
transmissão dos impulsos nervosos em todo o sistema nervoso e, se não houver 
uma pausa em determinados instantes, poder-se-á ter uma super transmissão de 
impulsos nervosos, gerando conseqüências drásticas como, por exemplo, batimento 
cardíaco e respiração não sincronizados até uma parada respiratória central, 
levando o indivíduo à morte. Felizmente, a acetilcolinesterase age hidrolisando a 
acetilcolina, regulando assim a transmissão dos impulsos nervosos no organismo 
vivo. Esta reação ocorre no interior do sítio ativo da acetilcolinesterase. 

1 .3.2.1 Interações da Acetilcolina com o Sítio Ativo da Acetilcolinesterase 
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Existem duas áreas importantes a serem consideradas: a região aniônica e a 
região catalítica (FIG. 1.4). 



FIG. 1 .4 - Região do sítio ativo da acetilcolinesterase, onde ocorrem as principais 
reações e interações com inibidores ou substratos (PATRICK, 2002; ORDENTLICH, 

1998). 

Pode-se notar que a acetilcolina interage com a acetilcolinesterase por uma 
ligação iônica com um glutamato (Glu) e uma ligação hidrogênio com uma tirosina 
(Tyr). Os resíduos histidina (His) e serina (Ser), no sítio catalítico, participam de um 
mecanismo de hidrólise. 

A região aniônica da acetilcolinesterase é muito similar à região aniônica dos 
receptores colinérgicos. Existe um poço hidrofóbico grande o suficiente para 
acomodar resíduos nas dimensões de radicais metil. A carga positiva do nitrogênio 
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pode interagir tanto com a carga negativa do Glu como também com aminoácidos 
aromáticos por interações dipolo induzidas. 


1 .3.2.2 Um dos Mecanismos de Hidrólise 


O resíduo histidina age como um catalisador ácido/base durante todo o 
mecanismo, enquanto a serina funciona como nucleófilo. Este não é o melhor papel 
da serina, tendo em vista que um álcool alifático é um nucleófilo fraco. De fato, a 
serina por si só é incapaz de hidrolisar um éster, porém, com a ajuda do resíduo 
histidina, isso se torna possível. O mecanismo se dá por várias etapas, como é 
mostrado na FIG. 1 .5. 
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FIG. 1.5 - Mecanismo de hidrólise da acetilcolina (PATRICK, 2002). 


Onde, na: 

Etapa 1 , a acetilcolina se aproxima, ligando-se a acetilcolinesterase. O resíduo 
serina age como um nucleófilo, usando um par de elétrons para formar um éster de 
acetilcolina. 
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Etapa 2, o resíduo histidina catalisa esta reação, agindo como uma base, 
removendo um próton. Isto torna a serina mais nucleofílica. 

Etapa 3, a histidina agora age como um catalisador ácido e protona o 
grupamento ‘OR’, tornando-o um melhor grupo de saída. 

Etapa 4, o grupamento carbonil expele o grupamento álcool do éster. 

Etapa 5, o grupamento acil da acetilcolina é covalentemente ligado ao sítio ativo. 
A colina deixa o sítio ativo, sendo substituída por água. 

Etapa 6, a água age como um nucleófilo e usa um par de elétrons do oxigênio 
para atacar o grupo acil. 

Etapa 7, como a água é normalmente um nucleófilo fraco, a histidina age 
novamente como um catalisador básico, capturando um próton. 

Etapa 8, a histidina agora age como um catalisador ácido, protonando o oxigênio 
da serina. 

Etapa 9, o grupamento carbonil é reconstituído e o resíduo serina é liberado. 
Etapa 10, o ácido etanóico deixa o sítio ativo e o ciclo recomeça. 

1.3.3 A Acetilcolinesterase Humana 

A partir de estudos de difração de raios-X, obteve-se a forma cristalográfica da 
acetilcolinesterase humana, complexada com fasciculina-ll (SOREQ.H.E, 1990). 

As estruturas de recombinantes de acetilcolinestarases selvagens e da mutante 
humana E202Q ( HuAChE ) complexadas com fasciculina-ll, uma toxina polipepitídica 
purificada que provém do veneno da cobra mamba verde, já foram relatadas. 
Posteriormente a estrutura do mutante E202Q foi refinada usando-se 2.7 A de 
resolução, obtendo-se assim a acetilcolinesterase a ser usada nesse trabalho, a 
qual se encontra depositada no servidor Protein Data Bank (PDB), sob o código 
1 B41 (KRYGER, 2000). 

1.3. 3.1 Inibição, Reativação e Envelhecimento da HuAChE 

1.3. 3. 1.1 Inibição 

A inibição da acetilcolinesterase ocorre no interior do sítio ativo, perto de um 
aminoácido serina específico, geralmente através da interação de um agente organo 
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fosforado, com o mesmo. No caso da acetilcolinesterase humana, o agente 
neurotóxico fosforilado liga-se covalentemente com o grupamento Ser-203. Esta 
ligação se dá entre o átomo de fósforo e o oxigênio da serina, com conseqüente 
desprendimento do grupo de saída do organofosforado (ORDENTLICH, 1996), 
como mostrado no mecanismo da FIG.1 .6. 
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FIG.1 .6 - Esquema geral de inibição da acetilcolinesterase, onde X é o grupo de 

saída do agente organofosforado. 


1.3. 3. 1.2 Reativação 


Para reativar a enzima inibida, geralmente se empregam substâncias conhecidas 
como oximas, FIG. 1.7, que competem fortemente com o aminoácido serina-203 
pelo organo fosforado, antes da etapa de envelhecimento. Uma das principais, e 
mais utilizada, tem sido a pralidoxima (EDDLESTON, 2002; KAMIL, 2004; TFIE 
MEDICAL LETTER, 2002). 
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FIG.1 .7 - Principais oximas empregadas atualmente. 

Essa etapa requer a hidrólise da ligação serina-fosforado, que pode ser feita com 
o auxílio do grupamento hidroxilamina, FIG.1 .8. 
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FIG.1 .8 - Hidrólise de organofosforados. 


A pralidoxima estando posicionada dentro do sítio ativo, FIG.1. 9, favorecerá a 
reação entre o grupo hidroxilamina e o éster fosforilado (PATRICK, 2002). 
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FIG.1. 9 - Reação entre o grupo hidroxilamina e o éster fosforilado. 

1.3. 3. 1.3 Envelhecimento 

A acetilcolinesterase inibida por neurotóxicos é considerada envelhecida quando 
há dealquilação da forma fosforilada. A enzima envelhecida dificilmente pode ser 
reativada por oximas. A janela terapêutica das oximas é, muitas das vezes, 
determinada pela taxa de envelhecimento. Na FIG.1. 10, a seguir, pode-se observar 
claramente todo o esquema de inibição, reativação e envelhecimento da 
acetilcolinesterase (EDDLESTON, 2002). 
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Com as facilidades computacionais disponíveis atualmente, tornou-se possível 
simular o comportamento fisiológico da AChE e compreender melhor as etapas de 
inibição, reativação e envelhecimento. Isso pode ser feito através do emprego das 
técnicas de modelagem e dinâmica molecular, entre outras. 
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FIG. 1.10- Inibição, reativação e envelhecimento da AChE. 


1.4 MODELAGEM MOLECULAR 

A modelagem molecular tem como principal objetivo reproduzir o comportamento 
de modelos de moléculas em sistemas moleculares. Até relativamente pouco tempo 
atrás, só se fazia modelagem de moléculas simples, usando modelos mecânicos ou 
lápis, papel e uma calculadora manual. O surgimento e a evolução das técnicas 
computacionais, todavia, revolucionaram a área e aumentaram sobremaneira o 
número e a complexidade dos sistemas que podem ser estudados por esta técnica 
(HEHRE, 1998). 

Há muita confusão sobre o significado dos termos: “química teórica”, “química 
computacional” e “modelagem molecular”. Certamente muitos aplicam os 03 termos 
para descrever aspectos de seus trabalhos de acordo com a ocasião. Como os 
modelistas moleculares utilizam todos esses métodos, não se deve preocupar com 
semântica e sim considerar qualquer técnica teórica ou computacional que forneça 
informações acerca do comportamento de sistemas moleculares como um exemplo 
de modelagem molecular (LEACH, 2001). 
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Quando se utiliza Modelagem Molecular, a primeira regra que se deve ter em 
mente é que nenhum método de cálculo é ideal para todas as aplicações, e que um 
grande esforço tem sido desenvolvido para achar métodos adequados às diferentes 
aplicações (HEHRE, 1998). 

Dentro da química quântica, os modelos Hartree-Fock (HF) são geralmente 
satisfatórios para uma extensa variedade de comparações termodinâmicas e 
cinéticas. Eles produzem um excelente resultado nas geometrias do equilíbrio e 
estado de transição, porém um razoável resultado para as conformações de 
equilíbrio. Existem limitações e ressalvas aos modelos Hartree-Fock. Primeiro, as 
interações eletrônicas instantâneas (correlações) não são consideradas devido à 
forma como a repulsão eletrônica é tratada. Isto significa que a posição instantânea 
de um elétron não é afetada pela presença de outro elétron em sua vizinhança. 
Também modelos HF não encontram muito sucesso na descrição de estruturas de 
compostos que contém metais de transição. Finalmente, deve ser levado em conta 
que os cálculos HF aumentam em custo computacional muito rapidamente com o 
aumento do tamanho da molécula. Isso faz com que eles não se tornem aplicáveis 
às moléculas com mais de 100 átomos. 

O insucesso dos modelos HF para descrever a geometria dos compostos que 
contém metais de transição pode ser conseqüência da incompleta descrição do 
acoplamento dos movimentos dos elétrons (“correlação eletrônica”). Modelos 
correlatos, que transpõem esta dificuldade, têm sido desenvolvidos. Dentre eles o 
mais popular é o MP2. Eles geralmente produzem excelentes descrições das 
geometrias do equilíbrio e estados de transição e conformações, além das 
grandezas cinéticas e termodinâmicas das reações, inclusive daquelas com o 
surgimento ou quebra de ligações. Esses modelos, todavia, são bem mais caros do 
ponto de vista computacional que os modelos HF, e seu uso se restringe 
usualmente a moléculas com menos de 20 átomos pesados. 

Uma alternativa aproximada para o problema da correlação eletrônica é 
encontrada nos conhecidos modelos de Densidade Funcional (DF). Ao invés de 
aproximar a solução multieletrônica por uma composição de soluções unieletrônicas, 
os modelos DF explicitam os efeitos multieltrônicos através de um termo de 
correlação baseado num problema multieletrônico “idealizado”. Comparativamente, 
o custo computacional envolvendo DF aumenta menos significativamente com o 
aumento da molécula do que se comparado a MP2 ou HF. Porém, existem 
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ressalvas no modelo. Diferente dos modelos anteriores cujas correlações podem ser 
melhoradas, o modelo DF só pode ser melhorado trazendo o problema 
multieletrônico “idealizado” mais próximo do “real”, porém esta etapa não é tão 
óbvia. Adicionalmente, os modelos DF não são tão simples de se trabalhar, pois 
requerem integrações numéricas introduzindo incertezas nos cálculos, fazendo com 
que eles devam ser cuidadosamente monitorados e controlados. 

Os métodos que utilizam a Mecânica Quântica e que são um pouco menos 
precisos são métodos semi-empíricos, conhecidos também como Métodos dos 
Orbitais Moleculares, e nesta lista se incluem os populares MNDO, AM1 e PM3. Os 
modelos semi-empíricos, que podem ser aplicados a moléculas compreendendo 
acima de 200 átomos, têm tido muito sucesso no cálculo da geometria de equilíbrio, 
incluindo as geometrias dos compostos que contém metais de transição e 
compostos organometálicos; e moderado sucesso para o cálculo das geometrias 
dos estados de transição. Eles, porém, não são satisfatórios para validações 
termodinâmicas e cinéticas ou para indicações conformacionais. 

Já o método da Mecânica Molecular pode facilmente ser aplicado a moléculas 
contendo 1000 átomos pesados ou mais em sua estrutura. Recentemente 
desenvolvido com o uso de campo de forças, produzem descrições quantitativas das 
estruturas e conformações das moléculas orgânicas. São exemplos o MMFF94 
(CLARK, 1989), SYBIL (HALGREN, 1996), GROMOS (VAN GUSTEREN, 1987), 
CHARMM (BROOKS, 1983) e o OPLS/AA (JORGENSEN, 1996). 

Quanto maior a complexidade de um sistema e quanto mais longo o tempo 
envolvido no processo em estudo, mais aproximações são requeridas. Neste 
momento a abordagem ab initio deverá ser aumentada ou substituída por uma 
parametrização semi-empírica ou empírica do modelo usado. Neste caso, a 
modelagem molecular está baseada inteiramente na análise por similaridade de 
dados estruturais e químicos já conhecidos. Muitas predições de estruturas de 
proteínas por homologia pertencem a esta categoria (VAN DER SPOEL, 2002). 

Várias abordagens, cada qual com suas limitações e sucessos, têm sido 
aplicadas para predizer estruturas de proteínas a partir de suas seqüências, dentre 
elas pode-se citar, por exemplo, a modelagem comparativa, conhecida também 
como modelagem por homologia. Algum sucesso vem sendo obtido com este 
método, mas ele apresenta algumas limitações como a sua dependência com a 
qualidade do alinhamento da seqüência alvo com a seqüência molde e a 
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necessidade da existência de uma boa similaridade seqüencial com o molde (acima 
de 35%) (HIGGINS, 2001). Com isso, pode-se passar para a próxima etapa que é a 
construção de modelos por homologia, que já se encontra quase que totalmente 
automatizada e envolve as seguintes etapas: identificação do molde e alinhamento 
de seqüências de aminoácidos, geração das coordenadas do modelo, otimização e 
validação do modelo. Contudo, sempre se pode fazer passo a passo. 

Então, Identifica-se o molde e o alinhamento de seqüências de aminoácidos para 
a determinação da estrutura da proteína-alvo (PEITSCH, 1995; PEITSCH, 1996). 
Neste ponto de vista, a busca do melhor molde pode ser feita dentro de uma família 
protéica. Assim por exemplo, se o objetivo é modelar uma HuAChE e só se conhece 
sua seqüência de aminoácidos, o caminho racional seria buscar um molde 
adequado dentro do grupo das acetilcolinesterases armazenadas no PDB 
(KRYGER, 2000). 

Uma vez definido(s) o(s) molde(s) potencial(is), parte-se para a etapa final que é 
a geração de um alinhamento de seqüências-alvo com as seqüências de todos os 
possíveis moldes. Se este alinhamento for feito entre proteínas de uma mesma 
família, crescem as chances de se obter maior precisão no modelo final. Esta 
envolve a modelagem das regiões estruturalmente conservadas, das alças e das 
cadeias laterais. Onde, no caso particular deste trabalho, isso foi feito de forma 
automatizada, através do seguinte procedimento: 

• Alinhamento da seqüência no formato FastA (PEARSON, 1988) com a 
proteína. Isso pode ser feito com o auxílio de programas como o Swiss- 
PdbViewer (GUEX, 2001). 

• O servidor SWISS-MODEL (PEITSCH, 1995; PEITSCH, 1996, GUEX, 1997), 
através de uma grande variabilidade de seqüência de alças contidas em seu 
banco de dados, completa os resíduos perdidos da proteína a qual deseja-se 
modelar. 

• É feita uma completa minimização de energia do conjunto, utilizando-se 
campo de forças apropriados. 

• Finalmente, o modelo completo é retornado como um arquivo com extensão 
pdb (Protein Data Bank). 

1.5 MECÂNICA MOLECULAR 
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Muitos dos sistemas com os quais se lida em modelagem molecular são 
demasiadamente grandes para serem abordados pela mecânica quântica ou por 
métodos semi-empíricos (HEHRE, 1998). Estes métodos levam em conta os 
movimentos dos elétrons de um sistema e mesmo que alguns sejam ignorados 
(como acontece nos métodos semi-empíricos) um grande número de partículas 
deve ainda ser considerado e os cálculos acabam consumindo um grande tempo 
computacional que os torna inviáveis para sistemas grandes. Os métodos de campo 
de forças, também conhecidos como Mecânica Molecular (MM), ignoram este fator e 
calculam a energia de um sistema como uma função apenas das posições 
nucleares. Isto torna a MM um médodo adequado para lidar com sistemas contendo 
um número grande de átomos. Em alguns casos, campos de força podem levar a 
respostas acuradas em um tempo computacional menor (HIGGINS, 2001). A MM, 
todavia, não pode prever propriedades que dependem da distribuição eletrônica em 
uma molécula tais como estados de transição ou distribuições de cargas. Para a MM 
funcionar corretamente é necessário fazer uma série de aproximações. A primeira 
delas é a aproximação de Born-Oppenheimer, sem a qual seria impossível 
descrever a energia como uma função de coordenadas nucleares. A MM é baseada 
num modelo mais simples das interações dentro de um sistema com contribuições 
de processos tais como o estiramento de ligações, o abrir e fechar de ângulos e as 
rotações em torno de ligações simples. Transferabilidade também é um atributo 
chave de um campo de forças, uma vez que ela permite que um conjunto de 
parâmetros desenvolvidos e testados em um número relativamente pequeno de 
casos seja aplicado a uma faixa muito mais ampla de problemas. Além disso, 
parâmetros desenvolvidos a partir de dados de pequenas moléculas podem ser 
usados para estudar moléculas muito maiores como as estruturas poliméricas, por 
exemplo. 

A maioria dos campos de força da MM em uso hoje em dia para sistemas 
moleculares, pode ser interpretada em termos de um simples conjunto de quatro 
componentes correspondentes às forças intra e intermoleculares dentro do sistema 
(VAN DER SPOEL, 2002). Penalidades energéticas estão associadas ao desvio de 
ângulos e ligações de seus valores de referência ou de equilíbrio; há uma função 
que descreve como a energia varia quando as ligações são giradas e, por fim, o 
campo de forças contém termos que descrevem a interação entre partes não 
ligadas do sistema. Campos de forças mais sofisticados podem incluir termos 
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adicionais, mas eles invariavelmente contêm estas quatro componentes. Uma 
característica atrativa desta representação é que os vários termos podem ser 
relacionados a variações em coordenadas internas específicas tais como 
comprimentos de ligações, rotação de ligações ou movimentos de átomos em 
relação a outros átomos. Isto torna mais fácil o entendimento de como variações 
nos parâmetros do campo de força afetam sua performance e também ajudam no 
processo de parametrização. Uma forma funcional deste campo de força é mostrada 
na EQ.1.1. 
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Em que: 

\>(r N ) é a energia potencial total, que é uma função das posições (r) de N partículas 
(normalmente átomos). O primeiro termo representa as interações entre pares de 
átomos ligados, representados aqui pelo potencial harmônico que dá o aumento na 
energia quando o seu comprimento n desvia do valor de referência n, 0 . O segundo 
termo é uma soma sobre todos os ângulos de valência (ângulos A-B-C) na 
molécula, novamente modelada usando um potencial harmônico. O terceiro termo é 
o potencial torcional, de caráter periódico, que modela como a energia varia quando 
as ligações giram. A quarta contribuição é o termo não ligado. Este é calculado entre 
todos os pares de átomos (i e j) que estão em diferentes moléculas ou na mesma 
molécula, mas separados por no mínimo 03 ligações (ou seja: tem uma relação 1,n 
onde n>4). Em um campo de força simples o termo não ligado é usualmente 
modelado usando um termo de potencial de Coulomb para interações eletrostáticas 
e um termo de potencial de Lennard -Jones para as interações de van der Waals. 


1.5.1 Termo para a Energia Potencial de Ligação 
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Envolvendo a interação entre dois átomos ligados por meio de um potencial 
harmônico (lei de Hooke), este termo reproduziria o movimento vibracional entre 
pares de átomos e tem a forma representada na EQ. 1 .2. 


v, 



( 1 . 2 ) 


Onde N b é o número total de pares de átomos ligados e para cada par são definidos 
os valores específicos em que k bn é a constante de força específica, n é a distância 
entre o par de átomos e n, 0 é o comprimento de referência da ligação. 



FIG.1 .1 1 - Energia potencial de ligação para um par de átomos. r n oscila com o 
tempo segundo um potencial harmônico em torno da posição de equilíbrio r on . 

1.5.2 Termo para a Energia Potencial Angular 

O desvio dos ângulos de seus valores de referência é também freqüentemente 
descrito, EQ.1.3, usando a lei de Hooke: 

, * < 13 > 

^ n = 1 

em que N e é o número total de trios de átomos ligados, k 0n a constante de força 
angular, 0 O o ângulo de equilíbrio e 0 a variação angular. 

A contribuição de cada ângulo é caracterizada por uma constante de força e um 
valor de referência. Como bem menos energia é necessária para distorcer um 
ângulo de sua posição de equilíbrio do que para estirar ou comprimir uma ligação, 
os valores numéricos das constantes de força são da ordem de 10 5 vezes menores. 
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FIG.1.12 - Energia potencial Angular para três átomos ligados. O ângulo entre eles, 
0 n varia harmonicamente em torno do ângulo de equilíbrio 0 on . 

1.5.3 Termo para o Potencial Torcional Próprio 

Termos para estiramento de ligação e deformação angular são geralmente 
conhecidos como termos de “graus de liberdade restritos” uma vez que quantidades 
substanciais de energia são requeridas para causar desvios significativos dos 
valores de referência. A maior parte das variações relativas na estrutura e nas 
energias é devido à inter-relações entre as contribuições torcionais e não ligadas. 
Muitos campos de forças são utilizados para modelar moléculas flexíveis nas quais 
as principais mudanças conformacionais são devidas a rotações em torno de 
ligações simples. Para simular isso, é essencial que o campo de forças represente 
apropriadamente os perfis de energia de tais mudanças. A maior parte dos campos 
de forças para moléculas orgânicas utiliza potenciais torcionais explícitos, também 
chamados de potenciais próprios, com uma contribuição de cada quarteto de 
átomos A-B-C-D ligados seqüencialmente no sistema (FIG.1.13). Os potenciais 
torcionais são quase sempre expressos como uma expansão em séries de co- 
senos. Uma forma funcional, EQ. 1.4, é: 


N 

u(4>) =^F^[ 1 + c °s(ot4> h -yj] 


n = 0 


(1.4) 


Onde (]) é o ângulo de torção; n é a multiplicidade ou o número de pontos mínimos 
na função quando a ligação é girada de 360°; yé o fator de fase que determina onde 
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o ângulo de torção passa por um mínimo e; V^ n é a constante que define a barreira 
de torção. 



FIG.1.13 - Energia potencial do torcional próprio. O ângulo O n varia em função da 
rotação do plano 1,2,3 em relação ao plano 2,3,4 ou vice versa. 

1.5.4 Termo para o Potencial Torcional Impróprio e Movimentos Angulares 
Fora do Plano 

Ou diedral impróprio, envolve quatro átomos não ligados sequêncialmente, 
sendo um central ligado aos outros três, sobre os quais se pode imaginar planos 
que compõem um ângulo diedral, que variaria harmonicamente. Este termo é 
fundamental para manter a planaridade dos átomos de hidrogênio, dos anéis 
aromáticos e da ligação pepitídica, assim como a quiralidade para certos átomos de 
carbono. 


v(u)=tk a M i -G ) o „) 2 

^ n = 1 


(1.5) 


Em que N u é o número total de quartetos de átomos envolvidos em diedrais 
impróprios, O a constante de força angular diedral impróprio, oo 0 n o ângulo diedral 
de equilíbrio e oo n sua variação angular. 
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FIG.1 .14 - Energia potencial do diedral impróprio. C0o n seria o ângulo entre os planos 
definidos pelos átomos 1 ,2,3 e os átomos 2,3,4. 

1.5.5 Termo para a Energia Potencial de Lennard-Jones 


Este termo inclui no campo de forças as interações de van der Waals entre 
átomos não ligados quimicamente. É um potencial de curto alcance (da ordem de Á) 
com um termo repulsivo, em virtude do princípio de exclusão de Pauli, que diz que 
dois elétrons não podem compartilhar os mesmos números quânticos, e um atrativo 
derivado das forças de dispersão, conhecidas também como forças de London, que 
surgem a partir de pequenos dipolos oscilantes criados pela flutuação da nuvem 
eletrônica. 
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Para um dado átomo / é computada sua energia de interação com os átomos j, 
de um total de Af, pois são excluídos seus primeiros e segundos vizinhos ligados, 
que são tratados pelos potenciais anteriores, ao passo que para os terceiros 
vizinhos parâmetros menores no potencial de Lennard_Jones são considerados. Os 
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parâmetros Ay e 8 V dependem dos tipos de pares de átomos (e se são terceiros 
vizinhos ou não) e nj é a distância entre eles. 



FIG.1.15 - Representação dos raios de van der Waals (R V dw), definidos como a 
metade da distância em que o potencial de Lennard-Jones do átomo vale zero 

(adaptada de WILTER, 2003). 


1.5.6 Termo para Energia Potencial Eletrostática 


Na definição dos parâmetros de carga para os átomos segundo um campo de 
forças, valores de cargas parciais lhe são atribuídos, de forma a se tentar reproduzir 
a diferença de eletronegatividade que eventualmente existe entre dois átomos 
ligados covalentemente. Pontes salinas e ligações hidrogênio, portanto, seriam 
decorrentes das interações eletrostáticas e teriam grande importância na 
estabilização das estruturas secundárias das proteínas e na especificidade entre 
ligante e proteína. 


t> 


el 


y My 

r tj 


(1.7) 


Na EQ. 1 .7, q , e g, são as cargas parciais para dois átomos quaisquer do sistema 
separados por uma distancia n/, e 0 é a permissividade do espaço livre, £ é a 
constante dielétrica do meio e, para um dado átomo /', AT é o número restante de 
átomos j, subtraindo seus primeiros e segundos vizinhos ligados. Contudo, o mais 
comum em DM é aplicar-se um raio de corte (Rc), isto é, dado um átomo i a partir 
do qual se delimita uma esfera de raio Rc, toma-se para cálculo de interação 
eletrostática e de van der Waals somente os átomos não ligados contidos nesta 
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esfera. Isto é feito fundamentalmente também para satisfazer as Condições 
Periódicas de Contorno (CPC), que se aplicam aos sistemas simulados. 



FIG.1.16 - Energia potencial de Coulomb entre dois átomos de cargas unitárias e 
opostas. Estão enfatizadas as linhas equipotenciais, sendo o tom da cor, indicador 
da intensidade do potencial, com o vermelho indicando valores positivos (átomo 
azul) e violeta para valores negativos (átomo vermelho) (adaptado de 

WILTER, 2003). 

Um sistema molecular seria mais precisamente descrito por uma abordagem que 
utilizasse a mecânica quântica. Contudo, face ao grande número de átomos 
envolvidos e de dificuldades computacionais para esta abordagem, as funções 
potenciais, de perfil empírico, são propostas de modelagem válidas derivadas de 
aproximações de teorias não clássicas e aplicam-se bem ao problema, uma vez que 
são calibradas a partir de informações experimentais e por cálculos quânticos sobre 
pequenas moléculas (BROOKS, 1983). 

Aceitando-se as limitações que um campo de forças clássico oferece, tais como, 
por exemplo, para o problema do congelamento “virtual” da água, quer seja pela 
ausência de potenciais mais apropriados para se tratar a energia de solvatação e 
outras interações intermoleculares, quer seja por questões de ordem entrópica não 
explicitamente tratadas, além da limitação que sistemas finitos impõem aos 
problemas envolvendo transição de fase, estes campos de forças respondem muito 
bem às simulações de dinâmica. 


1.6 OTIMIZAÇÃO DA ESTRUTURA 
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As forças interatômicas numa configuração inicial para uma estrutura molecular 
tridimensional, derivadas da EQ. 1.1, podem ser demasiadamente grandes, 
decorrendo grandes acelerações locais, implicando grandes velocidades e 
conseqüentes grandes deslocamentos que não satisfazem as condições para DM 
que, conseqüentemente, não funcionará. Assim, é necessária a minimização de sua 
energia, que é importante também para a remoção de qualquer tensão local. A 
EQ.1.1 gera uma hipersuperfície que apresenta mínimos locais, dos quais um ou 
mais podem ser globais (mínimos globais degenerados). Estes mínimos locais 
podem ser em números elevadíssimos já que dependem do número de graus de 
liberdade do sistema. E como este número é altíssimo para um sistema de 
macromoléculas biológicas, chegar-se a um mínimo global ao varrer-se 
completamente uma superfície multidimensional torna-se, hoje, praticamente 
inexeqüível. Portanto, o que se faz é procurar um mínimo ao se varrer uma parte da 
hipersuperfície potencial de forma a encontrar-se um ponto com a menor energia 
potencial local. O objetivo da consecução de um mínimo é ter uma conformação 
espacial que tenha relaxado as distorções nas ligações químicas, nos ângulos entre 
ligações e nos contatos de van der Waals. Um ponto mínimo (tanto local quanto 
global) é determinado quando as forças que atuam sobre o sistema são nulas. O 
mínimo global seria a configuração espacial para o sistema biológico cuja energia 
potencial total (EQ. 1.1) seria a menor dentre todas as outras, para qualquer ponto 
tomado na hipersuperfície. 

Apesar de existirem diversos métodos de otimização de geometria ou 
minimização de energia, abordar-se-á aqui apenas um dos mais largamente 
utilizado, o método de máximo declive ( steepest descent). 

1.6.1 Método do Máximo Declive 

O método do máximo declive (WIBERG, 1965) é um método de primeira 
derivada que converge vagarosamente nas proximidades do mínimo, mas é 
poderoso para configurações distantes de um mínimo de energia. De um modo 
geral, com este algoritmo é possível melhorar estruturas cristalográficas pouco 
refinadas, otimizar as construídas graficamente ou construir estruturas a partir de 
informações de RMN ou por homologia a estruturas conhecidas. 
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A força resultante sobre um átomo / é obtida pelo gradiente da energia potencial 
(dada pela EQ. 1 .8) em /: 


F,=-V£({r,}) 


( 1 . 8 ) 


O passo ou incremento nas coordenadas Ar, n de um átomo i, é dado na direção 
e sentido da força resultante sobre este átomo: 
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(1.9) 


F tn 

onde k n é um parâmetro de ajuste do tamanho do passo e é o vetor unitário na 

direção e sentido da força resultante sobre / no passo n. O algoritmo de minimização 
para o método do máximo declive é, então expresso na EQ.1.10, por: 


r i,n + 1 = r ,,n +K 
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( 1 . 10 ) 


Em que r vai dar a nova posição do átomo / no passo n+1 e k n é o parâmetro de 
ajuste do tamanho do passo. O que se obtém ao final da minimização é um conjunto 
de coordenadas { n } que, aplicado na EQ. 1.8, aproxima a força de zero, obtendo-se 
assim um mínimo local. A FIG.1.17 exemplifica o método para o caso bidimensional. 
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V(x,y) 


FIG.1 .17 - Ilustração do método de máximo declive para um dado campo de forças, 
cuja função (de duas variáveis apenas, para efeito didático) descreve a superfície 
que convergiu após 35 iterações. O caminho percorrido é uma geodésica, curva que 
cruza as linhas equipotenciais perpendicularmente ao longo dos gradientes 
máximos (adaptada de WILTER, 2003). 


1.7 SIMULAÇÕES POR DINÂMICA MOLECULAR 

Dinâmica molecular é uma técnica com a qual se estudam os movimentos em 
um sistema de partículas por simulação. Ela pode ser empregada tanto em sistemas 
de elétrons, átomos ou moléculas, como em sistemas macromoleculares (MARK, 
2000). Seus elementos essenciais são o conhecimento do potencial de interação 
entre as partículas e das equações de movimento que governam a dinâmica dessas 
partículas. O potencial pode variar do simples, como o gravitacional para interações 
entre estrelas, ao complexo, composto por vários termos, como o que descreve as 
interações entre átomos e moléculas. Para muitos sistemas, entre os quais os 
biomoleculares, as equações da dinâmica clássica são adequadas. Porém, para 
alguns problemas, como a evolução de galáxias, são incluídos efeitos relativísticos, 
enquanto que para outros, como reações químicas envolvendo tunelamento, são 
necessárias correções quânticas (KARPLUS, 1990; ISTVAN, 2001). 

O estado microscópico de um sistema pode ser especificado em termos das 
posições e momentos das suas partículas. Assim, pode-se escrever a Hamiltoniana 
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H de um sistema molecular clássico como a soma das energias cinética C e 
potencial V em função das séries de coordenadas generalizadas q , e de momentos 
generalizados p, de todos os N átomos do sistema (MARK, 2000): 


H({q„p l }) = C({p i }) + V({q i }) 


( 1 . 11 ) 


Onde q,=qi, q 2 ,..., qNat e p,=pi, p 2 , ..., p Nat 

A energia potencial V({qi}) contém os termos de interação inter e 
intramoleculares, de curto e longo alcance, e pode ser substituída pela função 
potencial V({ri}) da EQ. 1.1, tal que as coordenadas q-, sejam as coordenadas 
cartesianas n e p, seus momentos conjugados. A energia cinética assume a forma: 

( 1 . 12 ) 


Nat 

c({p,})=£ 


Onde m, é a massa do átomo /. 


A partir da H é possível construir as equações de movimento que governam a 
evolução temporal do sistema e suas propriedades dinâmicas. Como a energia 
potencial dada pela EQ. 1.1 independe das velocidades e do tempo, H é igual a 
energia total do sistema e as equações derivadas da Hamiltoniana: 



(1.13) 


dH 


Pi = - 


dq t 


(1.14) 


conduzem às equações do movimento de Newton: 
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(1.15) 



Pi = m t r t = 


MM 

dr f 


= Fi 


(1.16) 


em que, r, é a aceleração do átomo /, enquanto F é a força sobre o átomo / 
(GOLDSTEIN, 1980; BERENDSEN, 1986; ALLEN, 1987). 

A DM consiste, portanto, na resolução numérica das Equações de movimento de 
Newton para um sistema de interações de N átomos: 


ui 


av 

dt 2 




(1.17) 


em que, as forças são derivadas negativas duma função potencial v(n, r 2 ,...,r N ): 



As equações são resolvidas simultaneamente em pequenos passos de tempo. O 
sistema é seguido por algum tempo, mantendo-se a temperatura e pressão nos 
valores requeridos e as coordenadas são escritas em um arquivo de saída em 
intervalos regulares. As coordenadas como função de tempo representam a 
trajetória do sistema. Depois de mudanças iniciais, chega-se usualmente a um 
estado de equilíbrio dinâmico e, por uma média da trajetória de equilíbrio, muitas 
propriedades macroscópicas podem ser extraídas (VAN DER SPOEL, 2002). 

É útil, neste ponto, considerar as limitações da simulação da dinâmica molecular 
que reside, portanto, no fato de que para cada nanosegundo (ns) de simulação são 
necessários dois milhões de passos, dependendo do intervalo de tempo dado (dt) e, 
isso para uma macromolécula de 200 átomos poderia levar horas de processamento 
em um supercomputador, utilizando-se um algoritmo eficiente. Uma descrição e 
análise da eficiência de algoritmos para simulação de dinâmica molecular pode ser 
encontrada em Berendsen & van Gunsteren (BERENDSEN, 1986) e Allen & 
Tildesley (ALLEN, 1987). Este último inclui rotinas em linguagem FORTRAN para 
alguns métodos de simulação. 
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1.8 ANÁLISES ESTATÍSTICAS 


As simulações de DM geram uma quantidade enorme de dados cuja análise 
cuidadosa é fundamental para o sucesso do método. Os dados mais importantes 
gerados em arquivos são as posições, velocidades e acelerações de cada átomo do 
sistema ao longo do tempo. A partir deles outros resultados podem ser 
depreendidos. Contudo é prático, e muito relevante também, o registro das mais 
diversas energias presentes, tanto totais quanto de interação entre átomos e/ou 
grupos do sistema, além da temperatura, pressão, forças e outras propriedades 
físicas macroscópicas. As análises estatísticas utilizadas para esta dissertação 
foram à análise do desvio padrão e do gráfico de Ramachandran 
(RAMACHANDRAN, 1968). 

1.8.1 Desvio Padrão 

Ou Desvio da Raiz Média Quadrática (DRMQ) do inglês Root Mean Square 
Deviation (RMSD) (SPIEGEL, 1994), é o cálculo feito para comparação entre 
estruturas e, no caso dos resultados de DM, pode ser estendido à comparação entre 
as diversas conformações espaciais assumidas pelo sistema ao longo do tempo em 
relação à estrutura média temporal deste conjunto ou até de quadro a quadro 
quando se trata de um cálculo de DRMQ espacial, excluindo-se o solvente. Este 
parâmetro pode dar uma idéia do quanto a estrutura tridimensional flutuou ao longo 
do tempo, bem como permitir a observação de flutuações locais, como por exemplo, 
quais os resíduos com maior mobilidade durante uma dinâmica, para o caso de 
proteínas. O desvio padrão é dado pela E.Q.1.19, 


[T£ =7 (1 - 19) 

DesvioPadrão - V (X - X ) 2 
| n ;=1 

onde n é o número total de amostras, X, é a coordenada tridimensional de cada 
átomo da amostra i e x a posição média do átomo entre todas as amostras. Uma 
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boa fonte sobre desvio padrão e erro padrão em populações e amostragens pode 
ser vista em BOX et al . (BOX, 1 978) e Zar (ZAR, 1 999). 


1.8.2 Gráfico de Ramachandran 


Um importante indicador da qualidade estereoquímica de uma proteína é a 
distribuição dos ângulos torcionais <\> e 'P da cadeia principal. A distribuição destes 
pode ser examinada através do gráfico de Ramachandran (RAMACHANDRAN, 
1968). A partir da análise de dados cristalográficos de um conjunto enorme de 
proteínas, Ramachandran et al. (RAMACHANDRAN, 1968) mapearam regiões para 
os valores dos ângulos (j) e ‘P para resíduos de aminoácidos que seriam aceitáveis 
para a composição tridimensional da proteína, além de poder identificar regiões 
correspondentes referentes a estruturas secundárias, FIG.1.18. 


Legenda: 
a - Alfa Hélice. 

P - Folhas Beta. 

J - Regiões mais favoráveis. 
Q - Regiões menos favoráveis. 
| | - Regiões não favoráveis. 


” ISO 0 (j) (phi) 180° 



FIG.1.18 - Representação geral de um gráfico de Ramachandran 



1.9 OBJETIVOS DO PRESENTE TRABALHO 
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O objetivo deste trabalho, além de inicialmente constituir o modelo por homologia 
de duas alças incompletas da acetilcolinesterase humana, consistiu também em 
utilizar técnicas de modelagem e dinâmica molecular para estudar as interações da 
pralidoxima, assim como de seu análogo (deazapralidoxima) protonado e 
desprotonado dentro e fora do poço do sítio ativo da enzima, inibida pelo agente 
neurotóxico tabun (GA), buscando identificar as energias envolvidas nas 
simulações, além das trajetórias percorridas tanto pelo fármaco, quanto pelo 
análogo neutro e carregado negativamente, de forma que se possa sugerir melhores 
fármacos, tanto no caso de contaminação dos seres humanos por organofosforados 
empregados em guerra química como também por agrotóxicos, tendo em vista que 
são bastante similares. 


2 METODOLOGIA 
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2.1 MODELAGEM DA ACETILCOLINESTERASE HUMANA (HuAChE) E DOS 
COMPLEXOS ENVOLVIDOS NESTE TRABALHO 


2.1.1 Construção do Modelo 

A modelagem da HuAChE foi feita utilizando o servidor SWISS-MODEL 
(PEITSCH, 1995; PEITSCH, 1996; GUEX, 1997) e o programa Swiss-PdbViewer 
(GUEX, 1997). Inicialmente, pegou-se a estrutura tridimensional do cristal da 
HuAChE complexada com fasciculina-ll, depositada no servidor RCSB Protein Data 
Bank (BERMAN, 2000) pelo código 1B41 (KRYGER, 2000). Foi observado que a 
estrutura da HuAChE, já com as coordenadas da fasciculina-ll tirada manualmente 
do arquivo pdb, apresentava duas alças incompletas, a primeira entre os 
aminoácidos Pro492 e Pro495; a segunda entre os resíduos Pro258 e Asn265. 

Pegou-se a seqüência de aminoácidos em formato FASTA de um precursor da 
acetilcolinesterase humana, conhecido pelo código P22303, no servidor ExPASy 
(GASTEIGER, 2003). Foi feita a sobreposição entre as seqüências do P22303 e da 
1B41 sem a fasciculina-ll, através do programa Swiss-PdbViewer, de forma a tentar 
se obter um bom alinhamento seqüencial. Sendo assim, a melhor sobreposição foi 
entre os resíduos iniciais Asp36 da seqüência do precursor e Asp5 da HuAChE, 
finalizando-se com os resíduos Thr574 do precursor e Thr543 da HuAChE, 
respectivamente, como é mostrado na FIG.2.1. 
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FIG.2.1 Alinhamento entre os aminoácidos do precursor P22303 e da HuAChE. 


Após a etapa do alinhamento por sobreposição das seqüências em formato 
FASTA do precursor com o cristal, enviou-se o conjunto para o SWISS-MODEL, 
onde a partir de uma grande seqüência randômica de alças e aminoácidos, 
presentes no banco de dados do servidor, foram escolhidas as melhores para serem 
colocadas nas partes incompletas da enzima utilizando-se o programa ProModll 
(GUEX, 1997). Posteriormente foi feita uma minimização de energia, usando o 
campo de forças GROMOS 96 (VAN GUNSTEREN, 1996) e o método do máximo 
declive para, após retorno do modelo pronto, se fazer a validação por grafico de 
Ramachandran. 
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2.1.2 Determinação do Sítio Ativo da HuAChE 


A determinação dos resíduos do sítio ativo da HuAChE se deu a partir do 
alinhamento entre todos os aminoácidos desta com os aminoácidos da 
acetilcolinesterase do Torpedo Californica ( TcAChE ), uma espécie de peixe elétrico 
que possui abundante quantidade da enzima em seu órgão elétrico 
(SCHUMACHER, 1986), além de ter bastante similaridade com as 
acetilcolinesterases dos mamíferos e já possuir todos os resíduos de seu sítio ativo 
determinados (KOELLNER, 2000). Isso foi feito utilizando o software Swiss- 
PdbViewer (GUEX, 1997). A estrutura cristalográfica da TcAChE, encontra-se 
disponível no servidor RCSB Protein Data Bank (BERMAN, 2000), pelo código 
1EA5. 

2.1.3 Montagem e Parametrização dos Complexos HuAChE-GA-Reativador e 
HuAChE-GA - Análogos do Reativador 

Inicialmente, foram obtidas as coordenadas de menor energia da pralidoxima 
( 2PAM ) e do organofosforado tabun (GA), enantiômero P s , atracados dentro do sítio 
ativo da acetilcolinesterase do Torpedo Californica ( TcAChE ) conforme proposto por 
Castro et al. (CASTRO, 2002), utilizando-se o módulo Docking do programa Insightll 
2.3.5. 

Posteriormente, foi feita a sobreposição do modelo ( HuAChE ) com a TcAChE- 
GA-2PAM, utilizando o módulo Magic Fit, implementado no programa Swiss- 
PdbViewer (GUEX, 1997), tendo em vista que são bastante similares, de forma a se 
aproveitar tanto as coordenadas do GA e do 2PAM, salvando-se assim o conjunto 
das coordenadas em um único arquivo pdb, denominado HuAChE-GA-2PAM.pdb, 
observando-se que o GA foi ligado covalentemente ao aminoácido serina 203 da 
acetilcolinesterase humana, utilizando-se o mesmo programa. O mesmo foi feito 
para a deazapralidoxima (DZP) e para deazapralidoxima aniônica ( DZPanc ), onde a 
DZP é um análogo da 2PAM e lhe foi tirado o nitrogênio aromático, e a DZPanc é a 
mesma substância, só que agora, sem o hidrogênio ácido, de forma que ficou 
negativa. Estes foram colocados, através da sobreposição manual, exatamente em 
cima da pralidoxima. As estruturas dos compostos utilizados neste trabalho estão 
representadas na FIG.2.2. 
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FIG.2.2 - Fórmula estrutural plana dos compostos utilizados nesse trabalho, a) 
Organofosforado tabun (GA), b) Pralidoxima ( 2PAM ), c) Deazapralidoxima (DZP) e 

d) Deazapralidoxima aniônica ( DZP-anc ). 

Para implementar mais simulações neste trabalho, alterou-se manualmente as 
coordenadas das moléculas de 2PAM, DZP e DZPanc de forma que ficassem um 
pouco mais afastadas do ponto de partida, sendo assim, novas posições foram 
criadas, uma intermediária no interior do canal do sítio ativo e outra próxima à 
entrada do mesmo, dando um total de nove complexos. 

Tendo em vista que no campo de força GROMOS 87 (VAN GUNSTEREN, 
1987), implementado no programa GROMACS 3.2.1 (VAN DER SPOEL, 2002), não 
há parâmetros para compostos diferentes de aminoácidos, proteínas e ácidos 
nucleicos, houve a necessidade de parametrizar-se, independentemente, cada um 
dos compostos, inclusive o inibidor, seguindo o procedimento: 

• No caso da 2PAM, DZP e DZP-anc, as coordenadas (extensão pdb) foram 
enviadas para o servidor PRODRG (SCFIUETTELKOPF, 2004), onde a partir 
de um programa feito em Fortran, foram obtidos os arquivos de coordenadas 
(extensão gro) e de topologias (extensão itp), onde há parâmetros como, 
nome dos átomos, ângulos, ângulos de diedro, terceiros vizinhos e cargas 
atômicas, necessários para o reconhecimento do GROMACS 3.2.1 (VAN 
DER SPOEL, 2002) como novos resíduos. Isso tudo, de forma independente, 
para cada molécula dentro e fora do sítio ativo. 

• Já para o resíduo GA, além de se fazer todo o processo descrito 
anteriormente, teve-se que ligá-lo covalentemente ao resíduo serina 203 
(Ser203) da HuAChE. Essa etapa só foi possível após a inclusão de um novo 
resíduo, no arquivo aminoacids.dat (VAN DER SPOEL, 2002), onde estão os 
nomes de todos os resíduos presentes no banco de dados do GROMACS 
3.2.1 (VAN DER SPOEL, 2002), o qual se chamou de SGA (serina 203 ligada 
ao GA), além de incluir-se os parâmetros do arquivo de topologia (SGA. itp), 



58 



retomado pelo servidor PRODRG (SCHUETTELKOPF, 2004) no arquivo 
ffgmx.rtp (VAN DER SPOEL, 2002), onde estão as topologias de cada 
resíduo. Como nos arquivos de topologia retornados pelo servidor PRODRG 
(SCHUETTELKOPF, 2004), algumas cargas pontuais apresentavam-se 
nulas, fizeram-se otimizações de energia por mecânica quântica, utilizando o 
programa GAMESS US (SCHMIDT, 1993) na base Hartree Fock, 6-31 G* 
usando unidade de carga atômica CHELPG (BRENEMANN, 1990), para um 
melhor ajuste. 

• Foram preparados então, três conjuntos de três complexos diferentes para 
posteriores minimizações de energia e simulações de dinâmica molecular. 
São eles: HuAChE-GA-2PAM, HuAChE-GA-DZP e HuA ChE-GA-DZPanc, na 
posição proposta por Castro et al. (CASTRO, 2002), com uma distância 
aproximada de 3,00 A entre o fósforo do GA e o oxigênio da 2PAM, DZP e 
DZPanc e, mais dois conjuntos dos mesmos complexos, um com a distância 
entre os mesmos átomos, em torno de 5,00 A e o outro, 12,00 Â. 

O procedimento de parametrização de todas as moléculas usadas nesse 
trabalho está descrito com maior detalhe no apêndice A e, a adição de novas 
topologias no campo de forças GROMOS 87, no apêndice B. 

2.2 OTIMIZAÇÃO E DINÂMICA MOLECULAR DOS COMPLEXOS MONTADOS 

Para os complexos HuAChE-GA-(2PAM, DZP e DZPanc), descritos 
anteriormente, foram simuladas utilizando-se o pacote GROMACS 3.2.1 (VAN DER 
SPOEL, 2002), em caixas cúbicas de aproximadamente 722,855 nm 3 contendo 
21.332 moléculas de água e seis íons Na + para os complexos com a 2PAM, sete 
para os complexos com a DZP e cinco para os complexos com a DZPanc de forma 
a se neutralizar os sistemas. Estes sistemas foram minimizados utilizando-se o 
campo de forças GROMOS 87 (VAN GUNSTEREN, 1987), implementado no pacote 
GROMACS 3.2.1 (VAN DER SPOEL, 2002). O algoritmo de minimização 
empregado foi o do máximo declive (steepest descent) com critério de convergência 
de 418 kJ/mol/nm (10 kcal/mol/A). Após otimização, a HuAChE sem os ligantes 
(. 2PAM , DZP e DZPanc) foi validada através dos programas PROCHECK 
(LASKOWSKI, 1993) e WHAT IF (VRIEND, 1990) no site 
http://biotech.ebi. ac. uk:8400/. 
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Os complexos minimizados e validados foram, então, submetidos à dinâmica 
molecular em duas etapas: Inicialmente 50 ps de dinâmica a 300 K com restrição de 
todas as ligações da proteína para acomodação das moléculas de água no sistema, 
isso com condições periódicas de contorno. Após este procedimento, foram 
executadas nove simulações de dinâmica molecular de 1000 ps na temperatura de 
300 K, sem restrição de posições, tempo de integração de 2 fs, PME (DARDEN, 
1993; ESSMANN, 1995) para interações eletrostáticas, com raio de corte igual a 15 
A no espaço real, e um raio de corte de 15 A para interações de van der Waals. 

Os procedimentos para geração dos arquivos de entrada e execução da 
otimização e dinâmica de proteínas utilizando o pacote GROMACS 3.2.1 (VAN DER 
SPOEL, 2002) encontram-se descritos em detalhes no apêndice C. 

2.3 ANÁLISE DE ESTRUTURAS, GERAÇÃO DE FIGURAS E CÁLCULO DAS 
SUPERFÍCIES ELETROSTÁTICAS 

Para analisar as estruturas geradas após os cálculos de otimização e dinâmica 
molecular, foram utilizados os programas VMD (HUMPHREY, 1996) e PyMOL 
(WARREN, 2004). Os gráficos de variação de energia, da variação do DRMQ, das 
distancias entre átomos e da formação de ligações hidrogênio foram analisados com 
o programa XMGRACE (TURNER, 2004) e gerados no programa CALC 
(SINGLETON, 2005), sendo que este último foi baseado em critérios geométricos, 
utilizando-se um raio de corte variando de 0 a 0,35 nm e o ângulo formado entre o 
hidrogênio, doador e aceptor, variando de 0 a 60°. Os gráficos salsicha que 
descrevem o comportamento das enzimas após dinâmica foram gerados no 
programa MOLMOL (KORADI, 1996) e melhorados, para visualização, com o 
programa POVRAY (DAVID, 1999). 

Foram calculadas também, as superfícies eletrostáticas da 2PAM, DZP, DZPanc 
e do sítio ativo da HuAChE, a partir da solução numérica da equação de Poisson- 
Boltzmann (BAKER, 2000; BAKER, 2001; HOLST, 2000), um dos modelos mais 
populares para descrever interações eletrostáticas, utilizando-se o programa APBS 
(BAKER, 2004) que pode facilmente ser inserido como plugin no programa PyMOL 
(WARREN, 2004). 

2.4 RECURSOS COMPUTACIONAIS 
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O recurso computacional utilizado para executar os programas Swiss-PdbViewer 
(GUEX, 1997), PyMOL (WARREN, 2004), VMD (CADDIGAN, 2004), MOLMOL 
(KORADI, 1996), XMGRACE (TURNER, 1991), POVRAY (DAVID, 1999) e GAMESS 
US (SCHMIDT, 1993) foi um computador PC Duron 1,3 GHz com 640 MB de 
memória RAM e uma placa de vídeo NVIDIA RIVA TNT2 Model 64/ 32MB. Para a 
execução do pacote de dinâmica molecular GROMACS 3.2.1 foi utilizado um cluster 
de computadores composto por 01 PC Pentium IV de 2.4 GHz com 1GB de 
memória RAM e 04 PC Pentium IV de 2.6 GHz com 512 GB de memória RAM, 
disponível no laboratório de modelagem molecular do IME. 


3 RESULTADOS E DISCUSSÃO 

3.1 CONSTRUÇÃO DA HuAChE COMPLETA 
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A partir da sobreposição do cristal (código: 1B41) sem a fasciculina II, 
metodologia descrita na seção 2.1.1, obtém-se o modelo com as duas alças 
completas, como mostra a FIG. 3.1 . 



FIG.3.1 - Acetilcolinesterase humana com as duas alças completas. 

Quando se analisa este resultado, percebe-se com clareza que foi uma 
metodologia simples, tendo em vista que o alinhamento foi feito entre a sequência 
do cristal com a própria sequência completa da acetilcolinesterase humana, de 
forma que os aminoácidos ficassem bém alinhados. 

A etapa seguinte foi a validação da HuAChE, através do gráfico Ramachandran, 
para se conceber o quanto a modelagem foi boa, a partir do melhor valor dos 
ângulos Phi(O) x Psi(H J ), que descreve com clareza uma boa ou má localização dos 
aminoácidos em regiões denominadas de mais favoráveis, adicionalmente 
permitidas, generosamente permitidas e regiões não permitidas. A FIG. 3.2 mostra 
que se têm a maior parte dos resíduos da enzima localizados nas regiões mais 
favoráveis. 
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Phi ( graus ) 


Estatísticas do gráfico 


■ - Aminoácidos nas regiões mais favoráveis - 87,0 % 

O- Aminoácidos nas regiões adicionalmente permitidas- 1 1,7 % 

□ - Aminoácidos nas regiões generosamente permitidas- 1,4 % 

□ - Aminoácidos em regiões não permitidas - 0,0 % 

FIG.3.2 - Grafico de Ramachandran para todos os aminoácidos da HuAChE. 


Sendo assim, há 98,7 % dos resíduos em regiões mais permitidas, mostrando 
que este modelo é perfeitamente aceitável. 

Já os resíduos simbolizados por um pequeno triângulo e que não aparecem na 
estatística do gráfico, são meramente glicinas ou prolinas que, não têm importância 
em análises desse tipo, pelo fato de a glicina não possuir cadeia lateral, podendo 
assumir qualquer virtualmente quaisquer valores para os ângulos Phi(O) x Psi(’4*), e 
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a prolina, pelo contrário, possuindo uma cadeia lateral virtualmente rígida, 
delimitando fortemente a variabilidade dos referidos ângulos. 

3.2 OBTENÇÃO DO SÍTIO ATIVO DA HuAChE 


Como a TcAChE já tem todos os aminoácidos de seu sítio ativo determinados 
(KOELLNER, 2000) e apresenta grande similaridade com a HuAChE, através da 
sobreposição das duas e do alinhamento seqüencial entre elas (FRANÇA, 2005), 
obteve-se prontamente o resultado como mostrado na FIG. 3.3. 
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FIG.3.3 - Alinhamento sequencial entre a TcAChE e a HuAChE. Em vermelho estão 
os aminoácidos iguais entre as duas sequências. 


Através deste alinhamento, pode-se observar com clareza que a maior parte dos 
resíduos são idênticos, representados por asterisco (*). Os representados por ponto 
(.) são similares e, aqueles representados por espaço ( ), são diferentes. Sendo 
assim, aumentou-se a confiabilidade de se obter os resíduos do sítio ativo da 
HuAChE por uma simples comparação com os aminoácidos do sítio ativo da 
TcAChE, além de se confirmar essa veracidade com alguns artigos presentes na 
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literatura (SHAFFERMAN, 1996; ISTVAN, 2001; ZRINKA, 2003) que citam alguns 
dos aminoácidos do sítio ativo da acetilcolinesterase humana. Com isso, obteve-se 
os resíduos do sítio ativo da HuAChE como mostra a TAB. 3.1 . 


TcAChE HuAChE TcAChE HuAChE TcAChE HuAChE 


Gln-69 

Gln-71 

Leu- 127 

Leu- 130 

Phe-331 

Phe-338 

Tyr-70 

Tyr-72 

Tyr-130 

Tyr-133 

Tyr-334 

Tyr-341 

Asp-72 

Asp-74 

Glu-199 

Glu-202 

Gly-335 

Gly-342 

Gln-74 

Leu-76 

Ser-200 

Ser-203 

Trp-432 

Trp-439 

Ser-81 

Thr-83 

Ala-201 

Ala-204 

His-440 

His-447 

Trp-84 

Trp-86 

Trp-233 

Trp-236 

Gly-441 

Gly-448 

Asn-85 

Asn-87 

Trp-279 

Trp-286 

Tyr-442 

Tyr-449 

Gly-117 

Gly-120 

Leu-282 

Leu-289 


Glu-450 

Gly-118 

Gly-121 

Ser-286 

Ser-293 



Gly-119 

Gly-122 

Ile-287 

Val-294 



Phe-120 

Phe-123 

Phe-288 

Phe-295 



Tyr-121 

Tyr-124 

Arg-289 

Arg-296 



Ser- 122 

Ser- 125 

Phe-290 

Phe-297 



Gly-123 

Gly-126 

Glu-327 

Glu334 



Ser- 124 

Ala-127 

Phe-330 

Tyr-337 




TAB. 3.1 - Aminoácidos do sítio ativo da HuAChE obtidos por comparação com os 
resíduos do sítio ativo da TcAChE. Em negrito, estão os resíduos diferentes. 


Com esse resultado, foi calculada a identidade entre os resíduos das duas 
enzimas, a partir da EQ. 3.1 , obtida a partir de um simples cálculo por regra de três. 


r. (3- 1 ) 

/ = — xlOO 

r. 


onde n é o número de resíduos iguais, r t o número de resíduos totais de uma só 
seqüência e/a identidade que, neste caso, foi de aproximadamente 84%. 

3.3 CONSTRUÇÃO DOS COMPLEXOS 

3.3.1 Posição de Referência 

Inicialmente, foram usadas as coordenadas de um dos modelos proposto por 
Castro e colaboradores (CASTRO, 2002), que foi a TcAChE, ligada covalentemente 
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ao GA e com a pralidoxima ( 2PAM ) na conformação de menor energia (E-anti-anti), 
FIG. 3.4, docada no interior do sítio ativo, próxima ao organofosforado GA. 



- Carbono 

- Hidrogênio 
I - Oxigênio 
| - Nitrogênio 


FIG. 3.4 - 2PAM na comformação de menor energia (E-anti-anti). 


Posteriormente, fez-se a sobreposição do modelo ( HuAChE ) com a TcAChE 
complexada com GA e 2PAM, de forma a se aproveitar tanto as coordenadas da 
pralidoxima quanto as do organofosforado para finalmente ligar o GA 
covalentemente ao resíduo serina 203 da HuAChE que corresponde ao resíduo 
serina 200 da TcAChE, obtendo-se o primeiro complexo, a Acetilcolinesterase 
humana ligada ao GA e complexada com o 2PAM nas coordenadas obtidas por 
Castro et al. (CASTRO, 2002), o qual se chamou de HuAChE-GA-2PAM(A), FIG. 
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- Oxigênio 

- Nitrogênio 

- Fósforo 


- Nitrogênio ^ £ - Âmmoaciaos oàsicds 


Colocou-se então, tanto a DZP quanto a DZPanc, através de uma sobreposição 
manual, nas mesmas coordenadas da 2PAM de modo que as distâncias entre os 
oxigénios de ambas e o fósforo do GA também ficasse em torno de 3,0 A. Assim se 
obteve mais dois complexos, um chamado de HuAChE-GA-DZP(A) e o outro de 
HuAChE-GA-DZPanc(A), lembrando que (A) representará daqui em diante a 
posição de referência para a construção dos próximos complexos. 

3.3.2. Posição Intermediária no Canal do Sítio Ativo 

Neste caso, afastou-se manualmente a 2PAM do GA, de forma que a mesma 
ficasse próxima à saída do sítio ativo, com uma distância entre o oxigênio da 2PAM 
e o fósforo do organofosforado em torno de 5,0 A, como mostra a FIG. 3.6, para a 
seguir se montar novos complexos a partir destas novas coordenadas. 









- Nitrogênio ▼ | - Aminoácidos ácidos 

- Fósforo - Aminoácidos básicos 


FIG. 3.6 - Sítio ativo da HuAChE destacado, com a pralidoxima na posição 
intermediária com uma distância entre o oxigênio da 2PAM e o fósforo do GA em 

torno de 5,0 A. 
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Sobrepôs-se então a DZP e a DZPanc na 2PAM, de forma que ficassem nas 
mesmas coordenadas, tendo assim mais três novos complexos. São eles: 
HuA ChE-GA-2PAM(B) , HuAChE-GA-DZP(B), HuAChE-GA-DZPanc(B). 

3.3.3. Posição na Entrada do Sítio Ativo 

Finalmente foram montados os complexos, analogamente ao caso anterior, com 
a DZP e a DZPanc nas mesmas coordenadas da 2PAM, só que agora com essas 
três moléculas colocadas, manualmente, exatamente na entrada do sítio ativo e 
centralizadas na cavidade, FIG. 3.7B, com uma distância entre o oxigênio da 2PAM, 
DZP ou DZPanc e do fósforo do GA, em torno de 12,0 A, FIG. 3.7A. 



FIG. 3.7 - Novas coordenadas para montagem de mais três novos complexos, com 
a 2PAM na entrada do sítio ativo (representação em superfície de van der Waals) à 
uma distância de aproximadamente 12,00 Â entre o oxigênio da 2PAM e o fósforo 
do GA, sendo (A) vista lateral e (B) vista superior. 


Assim, analogamente ao que foi feito nos dois itens anteriores, criou-se mais três 
complexos. São eles: HuAChE-GA-2PAM(C), HuAChE-GA-DZP(C), HuAChE-GA- 
DZPanc(C). 


3.4 SIMULAÇÕES POR DINÂMICA MOLECULAR 


Após minimização de energia, os nove complexos foram submetidos a 
simulações por Dinâmica Molecular utilizando o campo de forças GROMOS 87. Este 
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procedimento teve o objetivo de colocar cada um deles e seus respectivos ligantes 
nas mesmas condições fisiológicas e de também servir como uma forma adicional 
de validação e refino das estruturas dos modelos uma vez que uma estrutura 
conservada no tempo durante uma simulação de dinâmica molecular é essencial 
para se considerar um modelo aceitável. 

Inicialmente, foi executada uma dinâmica molecular de 50 ps, a 300 K com 
restrição de posições para todo o sistema, exceto as moléculas de água, com o 
objetivo de garantir o equilíbrio das moléculas do solvente em torno dos resíduos da 
proteína. Após este procedimento, foi executada uma dinâmica molecular de 1000 
ps, sem restrição de posições, a 300 K, aplicando 2 fs de tempo de integração, PME 
(DARDEN, 1993; ESSMANN, 1995) para interações eletrostáticas, com raio de corte 
de 15 A no espaço real e um raio de corte de 15 Â para interações de van der 
Waals. No total, 1000 conformações foram obtidas durante esta simulação. 

Ao final das dinâmicas moleculares dispunha-se de um conjunto com várias 
conformações para os sistemas em estudo; os chamados quadros de dinâmica 
molecular (frames). A próxima etapa, então, seria a verificação da qualidade das 
dinâmicas, ou seja, determinar se a HuAChE apresentou boa estabilidade 
conformacional nestas nove simulações. Uma das formas de se obter valores 
quantitativos referentes à estabilidade de sistemas como os apresentados aqui é 
calculando-se o DRMQ (Desvio da Raiz Média Quadrática) para o esqueleto da 
estrutura da proteína em relação à estrutura média do conjunto. Os resultados que 
mostram as flutuações conformacionais estão ilustrados inicialmente aqui somente 
para o primeiro complexo, veja a FIG. 3.8, onde A representa o DRMQ em duas 
dimensões e B representa o DRMQ 3D. 

As avaliações de DRMQ na FIG. 3.8A foram feitas tomando-se para cada 
sistema 500 amostras, uma a cada 2 ps de simulação. Para a FIG. 3.8B, a 
amostragem nas simulações foi feita a cada 20 ps. Da FIG. 3.8A pode-se ainda 
observar que a média do DRMQ raramente ficou acima de 0,2 nm (2 A) o que pode 
ser considerado um valor bastante razoável quando comparado aos valores de 
resolução de cristalografia por raios X, cujos melhores resultados são da ordem de 
1,0 a 1,5 Â. Os picos vistos no gráfico indicam os pontos (resíduos) onde houve 
maior flexibilidade relativa para as vizinhanças destes referidos pontos. 
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FIG. 3.8 - (A) Gráfico da variação de DRMQ por resíduos para o complexo HuAChE- 
GA-2PAM(A). No eixo X está a seqüência de resíduos de aminoácidos para a 
cadeia da HuAChE, indo do resíduo 5 ao 544, onde os resíduos destacados na linha 
vermelha e na linha verde, foram obtidos através da modelagem por homologia. (B) 
Gráfico da variação de DRMQ 3D por resíduos para o complexo HuAChE-GA-2PAM 
(A), onde as setas vermelha e verde apontam as duas alças que possuem maiores 
variações do DRMQ, destacados respectivamente pelas mesmas cores no gráfico 

(A). 


Como a maioria dos picos máximos nunca ultrapassam 0,30 nm (3,0 Â), conclui- 
se que o modelo fica estabilizado durante as dinâmicas. O gráfico para a variação 
da energia total, apresentado na FIG. 3.9, corrobora estes resultados e mostra que 
a energia total para a proteína tende a estabilizar-se em torno dos 300 ps de 
dinâmica. 
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FIG. 3.9 - Gráfico de variação da energia total para o modelo complexo HuAChE- 
GA-2PAM(A) ao longo de 1000 ps de dinâmica molecular. 


Decidiu-se fazer então as análises dos gráficos de energia total para todos os 
outros complexos, a fim de se ter um critério de comparação, observando-se o 
comportamento da energia em todos eles ao longo das simulações. 

3.4.1 Análise dos Gráficos da Variação da Energia Total 

A idéia desta análise é verificar que a variação de energia total estabiliza-se, 
aproximadamente, a partir do mesmo período tempo para todas as simulações. 

Inicialmente, montou-se dois gráficos de energia total para, respectivamente, os 
complexos HuAChE-GA-DZP(A) e HuAChE-GA-DZPanc(A), que podem ser vistos 
na FIG. 3.10. Quanto ao gráfico da variação de energia total para o complexo 
HuAChE-GA-2PAM(A), este já foi mostrado anteriormente, FIG. 3.9. 

Posteriormente, fez-se o mesmo para os complexos HuAChE-GA-2PAM(B), 
HuAChE-GA-DZP(B) e HuAChE-GA-DZPanc(B), FIG. 3.11, finalizando-se com os 
complexos HuAChE-GA-2PAM(C), HuAChE-GA-DZP(C) e HuAChE-GA-DZPanc(C), 
gráficos estes, mostrados na FIG. 3.12. 
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Tempo ( ps ) 


FIG. 3.10 - Variação de energia total para os complexos (A) HuAChE-GA-DZP(A) e 

(B) HuAChE-GA-DZPanc(A). 
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FIG. 3.1 1 - Variação de energia total para os complexos (A) HuAChE-GA-2PAM(B), 
(B) HuAChE-GA-DZP(B) e (C) HuAChE-GA-DZPanc(B). 
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FIG. 3.12 - Variação de energia total para os complexos (A) HuAChE-GA-2PAM(C), 
(B) HuAChE-GA-DZP(C) e (C) HuAChE-GA-DZPanc(C). 

A partir das análises da FIG. 3.9 até a FIG. 3.12, pode-se verificar então que 
realmente todos os nove sistemas tendem a estabilizar-se, aproximadamente, a 
partir dos 300 ps de dinâmica molecular. 
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3.4.2 Análise dos Gráficos de Desvio Padrão 


Foram feitos cálculos de DRMQ temporal, em relação à média, sobre todos os 
átomos dos complexos para 1000 quadros, um a cada 1 ps, perfazendo assim, um 
total de 1 ns. Isso foi feito, separadamente para dois conjuntos em cada cálculo, um 
para a HuAChE+GA em relação a HuAChE+GA e o outro para a 2PAM, DZP ou 
DZPanc em relação à 2PAM, DZP ou DZPanc, dependendo do complexo. 
Considerando-se que o complexo poderia flutuar na caixa, cada quadro é ajustado 
pelo método dos mínimos quadrados a seu quadro precedente para efeito de 
cálculo de desvio padrão. Os resultados podem ser observados, desde a FIG. 3.13 
até a FIG. 3.15. 

Como se pode observar, em todos os três conjuntos de três complexos, para a 
proteína ligada ao organofosforado, o DRMQ tende a estabilizar-se logo nos 
primeiros pico segundos, não ultrapassando 0,8 Â. 

No caso da 2PAM pode-se notar a estabilização para os 3 complexos e para 
HuAChE-GA-2PAM(A), HuAChE-GA-2PAM(B) e HuAChE-GA-2PAM(C) pode-se 
notar também uma estabilização. Neste caso o DRMQ praticamente não passa de 
0,2 A. 

Quanto a DZP, a variação passa de 0,2 Â em alguns intervalos de tempo para o 
complexo HuAChE-GA-DZP(A)] passa de 0,4 Â para o complexo HuAChE-GA-DZP 
(B) em alguns intervalos de tempo e; passa de 0,5 A para o complexo HuAChE-GA- 
DZP(B) em uma quantidade considerável de intervalos de tempo. 

Já para a DZPanc, a variação do DRMQ em todos os casos, ultrapassa com 
maior freqüência 0,5 Â ao longo das dinâmicas, referentes aos complexos HuAChE- 
GA-DZPanc(A), HuAChE-GA-DZPanc(B) e HuAChE-GA-DZPanc(C), levando-se a 
supor que a DZPanc deve está se movimentando mais, ou seja, ela pode está 
sendo atraída ou repelida pela enzima com velocidade e aceleração maiores. 
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FIG. 3.13 - Gráfico de DRMQ temporal para as simulações: (A) HuAChE-GA-2PAM 
(A), (B) HuAChE-GA-DZP(A) e (C) HuAChE-GA-DZPanc(A), onde em azul é o 
DRMQ da enzima inibida em relação a própria enzima inibida e, em vermelho, da 
2PAM, DZP ou DZPanc em relação a eles mesmos, respectivamente em (A), (B) e 

(C). 
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FIG. 3.14 - Gráfico de DRMQ temporal para as simulações: (A) HuAChE-GA-2PAM 
(B), (B) HuAChE-GA-DZP(B) e (C) HuAChE-GA-DZPanc(B), onde em azul é o 
DRMQ da enzima inibida em relação a própria enzima inibida e, em vermelho, da 
2PAM, DZP ou DZPanc em relação a eles mesmos, respectivamente em (A), (B) e 

(C). 
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FIG. 3.15 - Gráfico de DRMQ temporal para as simulações: (A) HuAChE-GA-2PAM 
(C), (B) HuAChE-GA-DZP(C) e (C) HuAChE-GA-DZPanc(C), onde em azul é o 
DRMQ da enzima inibida em relação a própria enzima inibida e, em vermelho, da 
2PAM, DZP ou DZPanc em relação a eles mesmos, respectivamente em (A), (B) e 

(C). 
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O DRMQ 2D espacial, em relação à posição média, sobre cada resíduo de 
aminoácido da HuAChE foi calculado no intervalo de tempo de 1 a 1000 ps, de 1 em 
1 ps (1000 quadros) e, de 20 em 20 ps (50 quadros), para o DRMQ 3D, em relação 
ao primeiro quadro, para todas as nove simulações, foi calculado com o propósito 
de se observar quais regiões a HuAChE apresentam maiores flutuações e, se estas, 
são consideravelmente relevantes. Nas FIGS. 3.16 a 3.18, são apresentados esses 
gráficos. No caso do DRMQ 2D, tem-se em vermelho, a variação total da proteína 
em relação a ela mesma, em anis, somente a variação para da cadeia lateral e, em 
preto, somente a variação para o esqueleto peptídico (N, Ca, C). 

Para cada DRMQ 2D, em conjunto tem-se um DRMQ 3D, a fim de ter-se uma 
visão qualitativa de todas as regiões ao longo das dinâmicas. Maiores valores de 
DRMQ e maiores espessamentos dos tubos nas imagens indicam os resíduos que 
sofreram maiores variações. É comum notar que o DRMQ para os resíduos da 
região do sítio ativo e das alfa hélices e folhas beta, com valores mais baixos, 
revelam regiões mais estáveis. Poucos são os picos que excedem os 0,2 nm e mais 
raros ainda os que ultrapassam os 0,4 nm. Os picos que têm variações maiores são 
justamente pertencentes às regiões dos resíduos 259-264 e 493-494, alças já 
descritas na FIG. 3.8 anteriormente, ou então a qualquer outra alça localizada na 
periferia da enzima. 
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FIG. 3.16 - Gráficos do DRMQ espacial para os complexos: (A) HuAChE-GA-2PAM 
(A), (B) HuAChE-GA-DZP(A) e (C) HuAChE-GA-DZPanc(A). (A’), (B’) e (C’) são os 
cálculos de DRMQ 3D para os respectivos complexos em (A), (B) e (C). 
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FIG. 3.17 - Gráficos do DRMQ espacial para os complexos: (A) HuAChE-GA-2PAM 
(B), (B) HuAChE-GA-DZP(B) e (C) HuA ChE-GA-DZPanc(B ) . (A’), (B’) e (C’) são os 
cálculos de DRMQ 3D para os respectivos complexos em (A), (B) e (C). 
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FIG. 3.18 - Gráficos do DRMQ espacial para os complexos: (A) HuAChE-GA-2PAM 
(C), (B) HuAChE-GA-DZP(C) e (C) HuAChE-GA-DZPanc(C). (A’), (B’) e (C’) são os 
cálculos de DRMQ 3D para os respectivos complexos em (A), (B) e (C). 
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3.4.3 Análise das Possíveis interações da 2PAM, DZP E DZPanc com Resíduos 
do Sítio Ativo da HuAChE e o GA 

Após os cálculos de DRMQ espacial, decidiu-se então fazer análises das 
possíveis interações do fármaco ( 2PAM ) e de seu análogo protonado {DZP) e 
aniônico ( DZPanc ) com os aminoácidos do sítio ativo da HuAChE e com o 
organofosforado Tabun {GA). Sendo assim, a primeira análise feita foi da distância 
entre o fósforo do Tabun e os oxigénios da pralidoxima e da deazapralidoxima em 
suas formas neutra e aniônica. 

3.4.3. 1 Análise das Distâncias entre o Fósforo do GA e os Oxigénios da 2PAM, 
DZP e DZPanc 

Tendo em vista que a primeira interação que ocorre entre as oximas e a 
acetilcolinesterase inibida por agentes neurotóxicos fosforados é entre o oxigênio 
das oximas e o fósforo dos organofosforados, decidiu-se primeiramente, calcular as 
distâncias entre o fósforo do GA e os oxigénios da 2PAM, DZP e da DZPanc, para 
cada complexo. 

A partir da posição inicial proposta por Castro et al. (CASTRO, 2002), obteve-se 
o resultado mostrado na FIG. 3.19. Os gráficos referentes às coordenadas em que a 
distâncias iniciais foram, respectivamente, de 5 Â e de 12 A são mostrados nas FIG. 
3.20 e FIG. 3.21. 

Da FIG. 3.19, a impressão inicial é que a 2PAM estabiliza-se, aproximadamente, 
a partir dos 100 ps, em aproximadamente 5 A e a DZP em torno de 4 Â. Já a 
DZPanc, parece afastar-se, chegando aos 10 Â entre 200 e 500 ps, tendendo a 
aproximar-se novamente, chegando aos 7Â. 

Nas dinâmicas referentes à FIG. 3.20, tem-se a impressão destas distâncias 
ficarem estáveis, desde os 200 ps até os 800 ps. Posteriormente, a DZP afasta-se 
até aproximadamente 7 A, enquanto a DZPanc chega a 8 A, no final dos 1000 ps. Já 
a distância referente à 2PAM e o GA fica estabilizada em torno dos 6 Â. 

Para o caso da análise referente aos resultados apresentados na FIG. 3.21, a 
2PAM se aproxima nos primeiros 200 ps, estabilizando-se em torno de 8 Â ao longo 
da simulação, enquanto a DZP se afasta devagar, chegando aproximadamente aos 
15 A. Já a DZPanc se afastou bastante, chegando a 21 A. 
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Esses três conjuntos de três resultados, representados em cada gráfico, podem 
ser vistos também nas FIG. 3.22, FIG. 3.23 e FIG. 3.24, respectivamente, onde 
neste caso, fez-se as análises através de três quadros extraídos em instantes 
diferentes (0 ps, 500 ps e 1000 ps), com o propósito de ter-se também uma visão 
holística do ocorrido. Em cada uma dessas figuras, tem-se a representação dos 
resíduos do sítio ativo da HuAChE em superfície de van der Waals, com uma 
transparência de 70%, afim de se poder observar as estruturas da SGA, 2PAM, 
DZP e DZPanc , representadas em varetas, onde as distâncias são dadas em Â. 


2PAM DZP - DZPanc 



FIG. 3.19 - Distâncias entre o oxigênio da 2PAM e o fósforo do GA, em preto; 
distância entre o oxigênio da DZP e o fósforo do GA, em vermelho e; distância entre 
o oxigênio da DZPanc e o fósforo do GA, em azul. Cálculo referente aos complexos 
HuAChE-GA-2PAM(A), HuAChE-GA-DZP(A) e HuAChE-GA-DZPanc(A) 
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2PAM DZP - DZPanc 
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FIG. 3.20 - Distâncias entre o oxigênio da 2PAM e o fósforo do GA, em preto; 
distância entre o oxigênio da DZP e o fósforo do GA, em vermelho e; distância entre 
o oxigênio da DZPanc e o fósforo do GA, em azul. Cálculo referente aos complexos 
HuA ChE-GA-2PAM(B) , HuAChE-GA-DZP(B) e HuAChE-GA-DZPanc(B) 



FIG. 3.21 - Distâncias entre o oxigênio da 2PAM e o fósforo do GA, em preto; 
distância entre o oxigênio da DZP e o fósforo do GA, em vermelho e; distância entre 
o oxigênio da DZPanc e o fósforo do GA, em azul. Cálculo referente aos complexos 
HuA ChE-GA-2PAM( C) , HuAChE-GA-DZP(C) e HuAChE-GA-DZPanc(C) 
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FIG. 3.22 - Distâncias entre o fósforo do GA e os oxigénios da 2PAM, DZP e 
DZPanc, respectivamente (de cima para baixo), para as simulaões dos complexos 
HuAChE-GA-2PAM(A), HuAChE-GA-DZP(A) e HuAChE-GA-DZPanc(A) 
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FIG. 3.23 - Distâncias entre o fósforo do GA e os oxigénios da 2PAM, DZP e 
DZPanc, respectivamente (de cima para baixo), para as simulaões dos complexos 
HuAChE-GA-2PAM(B), HuAChE-GA-DZP(B) e HuAChE-GA-DZPanc(B) 
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FIG. 3.24 - Distâncias entre o fósforo do GA e os oxigénios da 2PAM, DZP e 
DZPanc, respectivamente (de cima para baixo), para as simulaões dos complexos 
HuAChE-GA-2PAM(C), HuAChE-GA-DZP(C) e HuAChE-GA-DZPanc(C) 


3.4. 3. 2 Análise dos Principais Aminoácidos em Posições Favoráveis às 
Interações com a 2PAM, DZP ou a DZPanc 


3.4.3.2.1 Nos complexos HuAChE-GA-2PAM(A), HuAChE-GA-DZP(A) e 
HuAChE-GA-DZPanc(A) 

No final da simulação referente ao primeiro complexo, observou-se que os 
aminoácidos Tyr341 e Tyr337 poderiam estar interagindo, ou em posições 
favoráveis a interações eletrostáticas com a 2PAM. Assim, foi decidido acompanhar 
a variação de algumas distâncias entre dois pares de átomos no decorrer da 
dinâmica. Os átomos são mostrados na FIG. 3.25A e a variação dessas distâncias 
na FIG. 3.25B, no decorrer de todo o tempo, onde ai é o oxigênio da Tyr341, a2 o 
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oxigênio da Tyr337 e a3, o nitrogênio positivo da 2PAM. Observou-se que as 
distâncias estabilizaram-se a partir dos primeiros 200 ps, uma em torno de 5 Â e a 
outra de 3 A 

Para o segundo complexo, o aminoácido mais próximo da DZP, no final da 
simulação, foi o Trp86, ficando numa posição extremamente favorável a interações 
do tipo empilhamento n (CAVALLI, 2004); interação que ocorre entre os orbitais n de 
dois ou mais anéis aromáticos. Calcularam-se então as distâncias entre 2 pares de 
carbonos pertencentes ao anel da DZP e ao anel do Trp86, chamados de a1-a2 e 
a3-a4. Os átomos em destaque e a posição propícia a interações desse tipo são 
exibidos na FIG. 3.26A e o cálculo dessas distâncias na FIG. 3.26B, onde se 
observou a estabilização destas em torno de 4 A, a partir dos 600 ps de simulação. 

Já para o terceiro complexo, observou-se que o oxigênio da DZPanc ficou 
suficientemente próximo ao hidrogênio da Tyr124, FIG. 3.27A, ao ponto de 
favorecer ligações hidrogênio (KYTE, 1995). Com isso, calculou-se essa distância 
ao longo da dinâmica, FIG. 3.27B, além de se calcular também o número de 
ligações hidrogênio entre a DZPanc e Tyr124, FIG. 3.27C. O resultado mostra que 
se forma pelo menos duas ligações hidrogênio durante quase todo o período de 
simulação, exceto no intervalo entre 350 e 500 os, aproximadamente, onde essa 
interação deixa de ocorrer devido a um afastamento entre o oxigênio da DZPanc e o 
hidrogênio da Tyr124. As duas ligações hidrogênio formam-se, provavelmente, 
devido a uma posição favorável entre o hidrogênio do grupamento hidroxila da 
Tyr124 e o nitrogênio da DZPanc. 
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FIG. 3.25 - (A) Último quadro da simulação referente ao complexo HuAChE-GA- 
2PAM(A), onde estão destacados os principais resíduos que podem interagir com a 
2PAM e as distâncias em A. (B) Distâncias calculadas ao longo da simulação. 
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FIG. 3.26 - (A) Último quadro da simulação referente ao complexo HuAChE-GA- 
DZP(A), onde estão destacados os principais resíduos que podem interagir com a 
DZP e as distâncias em Á. (B) Distâncias calculadas ao longo da simulação. 
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FIG. 3.27 - (A) Ultimo quadro da simulação referente ao complexo HuAChE-GA- 
DZPanc(A), onde estão destacados os principais resíduos que podem interagir com 
a DZPanc e a distância em Â. (B) Distância destacada em (A), calculada ao longo 
da simulação. (C) Número de ligações hidrogênio entre a DZPanc e a Tyr124. 
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3. 4. 3. 2. 2 Nos complexos HuAChE-GA-2PAM(B), HuA ChE-GA-DZP(B) e 
HuA ChE-GA-DZPanc(B) 


Para o primeiro complexo, ao final de 1000 ps de simulação, observou-se que os 
resíduos mais próximos e em condições favoráveis a algum tipo de interação com a 
2PAM foram o Trp86 e a Ser125, FIG. 3.28A, onde o anel aromático do triptofano 
ficou praticamente paralelo ao anel da 2PAM, podendo favorecer interações do tipo 
empilhamento n. Por outro lado, um dos oxigénios da Ser125 poderia está 
interagindo eletrostaticamente com o nitrogênio carregado positivamente da 2PAM. 
Decidiu-se então calcular essas distâncias de possíveis interações ao longo do 
tempo, FIG. 3.28B, onde se pode observar, uma estabilização dessas distâncias a 
partir dos primeiros 200 ps. 

No segundo complexo, o resíduo mais próximo da DZP foi uma Tyr337, FIG. 
3.29A, na ordem do comprimento de uma possível ligação hidrogênio. Calculou-se 
assim a variação desta ao longo do tempo, FIG. 3.29B, além de se computar se 
realmente estaria havendo interações desse tipo, FIG. 3.29C. O resultado nos 
mostra que houve uma efetiva estabilização a partir dos 800 ps e que realmente 
houve, ligações hidrogênio ao longo do tempo, exceto nos intervalos de, 
aproximadamente, 10 - 150 ps, 280 - 300 ps, 440 - 460 ps e 720 - 780 ps, 
corroborando com o resultado exibido na FIG. 3.29B. 

Já para o terceiro complexo, o grupamento hidroxila da DZPanc ficou próximo o 
suficiente da hidroxila da Tyr124 ao ponto de proporcionar também possíveis 
ligações hidrogênio, FIG. 3.30A. A distância praticamente não ultrapassou 2 A 
durante quase toda a simulação, FIG. 3.30B, favorecendo formação de duas 
ligações hidrogênio estáveis, durante todo o tempo, FIG. 3.30C. 
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FIG. 3.28 - (A) Último quadro da simulação referente ao complexo HuAChE-GA- 
2PAM(B), onde estão destacados os principais resíduos que podem interagir com a 
2PAM e as distâncias em A. (B) Distâncias calculadas ao longo da simulação. 
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FIG. 3.29 - (A) Último quadro da simulação referente ao complexo HuAChE-GA- 
DZP(B), onde estão destacados os principais resíduos que podem interagir com a 
DZPe a distâncias em Â. (B) Distância destacada em (A), calculada ao longo da 
simulação. (C) Número de ligações hidrogênio entre a DZP e a Tyr337. 
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FIG. 3.30 - (A) Último quadro da simulação referente ao complexo HuAChE-GA- 
DZPanc(B), onde estão destacados os principais resíduos que podem interagir com 
a DZPanc e a distâncias em A. (B) Distância destacada em (A), calculada ao longo 
da simulação. (C) Número de ligações hidrogênio entre a DZPanc e a Tyr124. 
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3. 4. 3. 2. 3 Nos complexos HuAChE-GA-2PAM(C), HuAChE-GA-DZP(C) e 
HuA ChE-GA-DZPanc(C) 


Para o complexo HuAChE-GA-2PAM(C), observou-se que ao final de 1000 ps de 
simulação, a pralidoxima ( 2PAM ) entrou no poço do sítio ativo da HuAChE, ficando 
ancorada em coordenadas favoráveis a interações do tipo empilhamento n com o 
Trp286 e possíveis interações eletrostáticas, respectivamente, com o Asp74 e com a 
Tyr341. Os átomos que interagem estão destacados na FIG. 3.31 A em preto, as 
distâncias calculadas ao longo da simulação na FIG. 3.31 B e a entrada da 2PAM no 
sítio ativo pode ser revista na FIG. 3. 24. 

O resultado mostra que, para possíveis interações pi-stacking a1-a2 e a3-a4, 
houve estabilização dessas distâncias após, aproximadamente, os primeiros 50 ps. 
Para a distância referente à interação eletrostática entre o oxigênio do grupamento 
OFI da Tyr341 e o nitrogênio neutro da 2PAM (a5-a6), começou a se estabilizar a 
partir dos 300 ps. Já para a distância entre o nitrogênio positivo da 2PAM e o 
oxigênio do grupamento carbonila do Asp74 (a7-a8), estabilizou-se a partir dos 650 
ps. 

No complexo HuAChE-GA-DZP(C), FIG. 3.32, observa-se que a DZP, quando 
colocada na entrada do sítio ativo, foi levemente repelida, ficando suficientemente 
distante dos resíduos do sítio ativo situados na periferia do canal, de forma a não 
favorecer mais nenhum tipo de interação com estes aminoácidos, no final dos 1000 
ps de simulação. 

Já no complexo HuAChE-GA-DZPanc(C), ao final dos 1000 ps, parece que a 
DZPanc foi mais fortemente repelida que a DZP no caso anterior, ficando mais 
distante ainda de possíveis interações com resíduos da entrada do canal, FIG. 
3.33A. Pode-se observar também que no final da simulação, o canal do sítio ativo 
ficou totalmente fechado, quando se analisou o sítio ativo em termos de superfície, 
FIG. 3.33B. 
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FIG. 3.31 - (A) Último quadro da simulação referente ao complexo HuAChE-GA- 
2PAM(C), onde estão destacados os principais resíduos que podem interagir com a 
2PAM e a distâncias em A. (B) Distâncias destacadas em (A), calculadas ao longo 

da simulação. 
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FIG. 3.32 - Vista superior do ultimo quadro referente ao complexo HuAChE-GA-DZP 
(C), onde estão destacados a DZP e os resíduos mais próximos dela na saída do 

sítio ativo com o SGA no fundo do poço. 
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- Oxigênio 

- Nitrogênio 
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FIG. 3.33 - Último quadro referente ao complexo HuAChE-GA-DZPanc(C) com o 
sítio ativo representado em superfície de van der Waals , a DZPanc e o SGA em 
varetas. Onde em (A) tem-se a vista lateral e em (B), a vista superior. 


3.4.3. 3 Cálculo das Superfícies Eletrostáticas 


Com o propósito de se tentar saber o “porquê” do afastamento tanto da DZP 
quanto da DZPanc e da aproximação da 2PAM pelo sítio ativo da enzima, quando 
colocados na entrada do canal, decidiu-se calcular a superfície eletrostática da 
2PAM, DZP, DZPanc e do sítio ativo da HuAChE. 
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A escala de energia potencial eletrostática é diferente para cada superfície e sua 
dimensão está em kcal/mol. Foram calculadas usando-se constante dielétrica igual a 
20, para os fármacos e para o sítio ativo, temperatura de 310 K e o potencial 
eletrostático com uma variação entre -9 e +4 kcal/mol para o sítio ativo da HuAChE, 
-4 e +4 kcal/mol tanto para a 2PAM quanto a DZPanc e, -1 e +1 kcal/mol para a 
DZP. Os resultados desses cálculos são mostrados na FIG. 3.34. 

O resultado mostra que praticamente toda a superfície da 2PAM é positiva, a 
superfície da DZP é predominantemente neutra, com poucas regiões positivas e 
negativas enquanto que, tanto a DZPanc quanto o sítio ativo da HuAChE, são 
predominantemente negativos. 



FIG. 3.34 - Cálculo da variação de energia potencial eletrostática para (A) 
2PAM, (B) DZP, (C) DZPanc e (D) sítio ativo da HuAChE. 
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4 CONCLUSÕES 


Na primeira etapa deste trabalho, foi feita a modelagem por similaridade de duas 
alças incompletas do cristal da acetilcolinesterase humana com a sua própria 
seqüência. O resultado se mostrou satisfatório, através de uma análise pelo gráfico 
de Ramachandran, para se poder partir para a etapa de simulação por dinâmica 
molecular. 

Ao se simular um meio biológico utilizando o modelo da acetilcolinesterase 
humana ( HuAChE ) utilizada neste trabalho junto com o fármaco e inibida pelo 
organofosforado GA, corroborou com o que já vem sendo descrito e proposto pela 
literatura, em termos das oximas serem os melhores fármacos. 

Os resultados obtidos neste trabalho mostram que: 

• Quando a 2PAM é colocada no interior do sítio ativo, tende a permanecer em 
uma posição favorável à interação com o GA. 

• Quando a DZP e a DZPanc são colocadas tanto nas coordenadas propostas 
por Castro et al. (CASTRO, 2002) quanto nas coordenadas intermediárias, no 
canal, tendem a permanecer dentro do sítio, possivelmente devido a 
interações eletrostáticas e interações do tipo empilhamento n com alguns 
aminoácidos. 

• Quando a 2PAM é colocada próxima à entrada do canal do sítio ativo, tende 
a entrar pelo fato de sua distribuição de cargas ser predominantemente 
positiva e a do sítio ativo ser predominantemente negativa. Já a DZP tende a 
sair e, a DZPanc, é repelida mais fortemente, devido a sua distribuição de 
cargas ser parecida com a do sítio ativo. 

Como continuação deste trabalho, são sugeridos os seguintes estudos: 

• Estudo por dinâmica molecular da reativação da acetilcolinesterase inibida 
por outros organofosforados, utilizando como fármacos, outras oximas, 
preferencialmente aquelas com mais de uma carga positiva. 

• Estudo por dinâmica molecular de oximas fosforiladas, dentro do canal do 
sítio ativo da acetilcolinesterase humana. 
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• Estudo da reativação da acetilcolinesterase por métodos híbridos de 
mecânica molecular/mecânica quântica (QM/MM), de forma a confirmar as 
hipóteses propostas neste trabalho. 

• Utilização do modelo completo da acetilcolinesterase humana, assim como as 
diversas informações obtidas nesta pesquisa, para projetar novos potenciais 
desinibidores. 

• Cálculos de energia livre e docking, para os fármacos estudados neste 
trabalho, em diferentes posições. 
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6.1 APÊNDICE A: PARAMETRIZAÇÃO DAS MOLÉCULAS 2PAM, DZP E DZPanc 


Parametrização, neste trabalho, é o preparo do arquivo com topologias para 
moléculas não contempladas pelos campos de forças ora presentes do pacote no 
GROMACS. 

O procedimento segue o roteiro: 

1- Acessa-se a página da internet: 

http://davapd .bioch.dundee.ac.uk/programs/prodrg/prodrg.html. 

2- Copia-se as coordenadas de todos os átomos da molécula, 
presentes no arquivo pdb, ao qual queremos parametrizar e cola-se 
na caixa de diálogos do servidor . 

3- Escolhe-se Yes para Chirality (quiralidade) e Full charge (carga total) 
e, No para Energy minimization (minimização de energia). 

4- Executa-se o programa através do botão RunPRODRG e espera-se 
por algum momento para que os arquivos com extensão itp e gro, 
respectivamente, de topologias e coordenadas, estejam disponíveis, 
em um arquivo compactado, para ser baixado. 

5- Faz-se o refinamento das cargas atômicas, presentes nos arquivos 
itp através de cálculos quânticos; neste trabalho usamos o GAMESS 
US. 

6- Copia-se os arquivos gro e itp para a pasta onde serão executados 
os programas do pacote GROMACS 3.2.1 . 

O conteúdo dos arquivos de topologia para a 2PAM, DZP e DZPanc são 
mostrados a seguir: 

2 P AM. itp : 

[ moleculetype ] 

; name nrexcl 
2PAM 3 


This file was generated by PRODRG version 040430.0567 
PRODRG written/copyrighted by Daan van Aalten 

Questions/comments to dava@davapcl.bioch.dundee.ac.uk 

When using this software in a publication, cite: 

D.M.F. van Aalten, R. Bywater, J.B.C. Findlay, 

M. Hendlich, R.W.W. Hooft and G. Vriend, 

PRODRG, a program for generating molecular topologies 
and unique molecular descriptors from coordinates 
of small molecules . 

Journal of Computer Aided Molecular Design (1996), 10, 255-262. 
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[ atoms ] 


nr 

type 

resnr 

resid 

atom 

cgnr < 

charge 

1 

CH3 

1 

2PAM 

C2 

1 

0.2927 

2 

NR6* 

1 

2PAM 

N 

1 

-0.4126 

3 

CR61 

1 

2PAM 

Cl 

1 

0.3683 

4 

CR61 

1 

2PAM 

C3 

1 

-0.0887 

5 

CR61 

1 

2PAM 

C6 

1 

0.3257 

6 

CR61 

1 

2PAM 

C4 

2 

-0.1690 

7 

CB 

1 

2PAM 

C 

2 

0.4666 

8 

CR61 

1 

2PAM 

C5 

2 

0.3547 

9 

N 

1 

2PAM 

NI 

2 

-0.2124 

10 

OA 

1 

2PAM 

O 

2 

-0.3879 

11 

HO 

1 

2PAM 

HAA 

3 

0.4626 


[ bonds ] 


ai 

aj 

fu 

cO 

cl 

1 

2 

1 

0 . 148 

334720 . 0 

2 

3 

1 

0 . 140 

334720 . 0 

2 

7 

1 

0 . 140 

334720 . 0 

3 

4 

1 

0.139 

418400 . 0 

4 

5 

1 

0.139 

418400 . 0 

5 

6 

1 

0.139 

418400 . 0 

6 

7 

1 

0.139 

418400 . 0 

7 

8 

1 

0.139 

418400 . 0 

8 

9 

1 

0.132 

418400 . 0 

9 

10 

1 

0.136 

376560.0 

10 

11 

1 

0.100 

313800.0 

pairs 

] 



ai 

aj 

fu 

cO 

cl 

1 

4 

1 



1 

6 

1 



1 

8 

1 



2 

5 

1 



2 

9 

1 



3 

6 

1 



3 

8 

1 



4 

7 

1 



5 

8 

1 



6 

9 

1 



7 

10 

1 



8 

11 

1 




0 . 148 

334720.0 ; 

C2 

N 

0 . 140 

334720.0 ; 

N 

Cl 

0 . 140 

334720.0 ; 

N 

C 

0.139 

418400.0 ; 

Cl 

C3 

0.139 

418400.0 ; 

C3 

C6 

0.139 

418400.0 ; 

C6 

C4 

0.139 

418400.0 ; 

C4 

C 

0.139 

418400.0 ; 

C 

C5 

0.132 

418400.0 ; 

C5 

NI 

0.136 

376560.0 ; 

NI 

O 

0.100 

313800.0 ; 

O 

HAA 


; 

C2 

C3 

; 

C2 

C4 

; 

C2 

C5 

; 

N 

C6 

; 

N 

NI 

; 

Cl 

C4 

; 

Cl 

C5 

; 

C3 

C 

; 

C6 

C5 

; 

C4 

NI 

; 

C 

O 

; 

C5 

HAA 


[ angles ] 


ai 

aj 

ak 

fu 

cO 

cl 






1 

2 

3 

1 

120.0 

418.4 

120.0 

418.4 ; 

C2 

N 

Cl 

1 

2 

7 

1 

120.0 

418.4 

120.0 

418.4 ; 

C2 

N 

c 

3 

2 

7 

1 

120.0 

418.4 

120.0 

418.4 ; 

Cl 

N 

c 

2 

3 

4 

1 

120.0 

418.4 

120.0 

418.4 ; 

N 

Cl 

C3 

3 

4 

5 

1 

120.0 

418.4 

120.0 

418.4 ; 

Cl 

C3 

C6 

4 

5 

6 

1 

120.0 

418.4 

120.0 

418.4 ; 

C3 

C6 

C4 

5 

6 

7 

1 

120.0 

418.4 

120.0 

418.4 ; 

C6 

C4 

C 

2 

7 

6 

1 

120.0 

418.4 

120.0 

418.4 ; 

N 

C 

C4 

2 

7 

8 

1 

120.0 

418.4 

120.0 

418.4 ; 

N 

C 

C5 

6 

7 

8 

1 

120.0 

418.4 

120.0 

418.4 ; 

C4 

C 

C5 

7 

8 

9 

1 

120.0 

418.4 

120.0 

418.4 ; 

C 

C5 

NI 

8 

9 

10 

1 

120.0 

418.4 

120.0 

418.4 ; 

C5 

NI 

O 

9 

10 

11 

1 

109.5 

397.5 

109.5 

397.5 ; 

NI 

O 

HAA 


[ dihedrals ] 

;ai aj ak al fu 

2 7 3 1 2 

7 2 8 6 2 


cO cl m cO cl m 
0.0 1673.6 0 0.0 1673.6 0 ; IDI 

0.0 1673.6 0 0.0 1673.6 0 ; IDI 


N C Cl C2 
C N C5 C4 
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2 

3 

4 

5 

2 

0. 

.0 

1673 

. 6 

0 

0 

. 0 

1673. 

. 6 

0 

IDI 

N 

Cl 

C3 

C6 

3 

4 

5 

6 

2 

0. 

.0 

1673 

. 6 

0 

0 

. 0 

1673. 

. 6 

0 

IDI 

Cl 

C3 

C6 

C4 

4 

5 

6 

7 

2 

0. 

.0 

1673 

. 6 

0 

0 

. 0 

1673. 

. 6 

0 

IDI 

C3 

C6 

C4 

C 

5 

6 

7 

2 

2 

0. 

.0 

1673 

. 6 

0 

0 

. 0 

1673. 

. 6 

0 

IDI 

C6 

C4 

C 

N 

6 

7 

2 

3 

2 

0. 

.0 

1673 

. 6 

0 

0 

. 0 

1673. 

. 6 

0 

IDI 

C4 

C 

N 

Cl 

7 

2 

3 

4 

2 

0. 

.0 

1673 

. 6 

0 

0 

. 0 

1673. 

. 6 

0 

IDI 

C 

N 

Cl 

C3 

2 

7 

8 

9 

1 

180. 

.0 

41 

. 8 

2 

180 

. 0 

41 . 

.8 

2 

Dl 

N 

C 

C5 

NI 

7 

8 

9 

10 

1 

180. 

.0 

33 

.5 

2 

180 

. 0 

33. 

.5 

2 

Dl 

C 

C5 

NI 

O 

8 

9 

10 

11 

1 

180. 

.0 

16 

. 7 

2 

180 

. 0 

16. 

. 7 

2 

Dl 

C5 

NI 

0 

HAA 


DZP.itp : 

[ moleculetype ] 
; name nrexcl 
DZP 3 


This file was generated by PRODRG version 040604.0507 
PRODRG written/copyrighted by Daan van Aalten 

Questions/comments to dava@davapcl.bioch.dundee.ac.uk 

When using this software in a publication, cite: 

A. W. Schuettelkopf and D. M. F. van Aalten (2004) . 
PRODRG - a tool for high-throughput crystallography 
of protein-ligand complexes. 

Acta Crystallogr. D60, in press. 


[ atoms ] 


nr 

type 

resnr 

resid atom 

cgnr 

charge 




1 

CH3 

1 

DZP CAJ 

1 

0 . 0299 




2 

CB 

1 

DZP CAF 

1 

0 . 0117 




3 

CR61 

1 

DZP CAE 

1 

-0 . 0527 




4 

CR61 

1 

DZP CAD 

1 

0 . 0101 




5 

CR61 

1 

DZP CAC 

2 

-0 . 0134 




6 

CR61 

1 

DZP CAB 

2 

0 . 0207 




7 

CB 

1 

DZP CAA 

2 

-0 . 004404 




8 

CR61 

1 

DZP CAG 

2 

0.3362 




9 

NL 

1 

DZP NAH 

2 

-0.2845 




10 

OA 

1 

DZP OAI 

3 

-0.5159 




11 

HO 

1 

DZP HAA 

3 

0 .4623 



bonds ] 







ai 

aj fu 

c0 

cl 





1 


2 1 

0.153 

334720 . 0 

0.153 

334720.0 ; 

CAJ 

CAF 

2 


3 1 

0.139 

418400 . 0 

0.139 

418400.0 ; 

CAF 

CAE 

2 


7 1 

0.139 

418400 . 0 

0.139 

418400.0 ; 

CAF 

CAA 

3 


4 1 

0.139 

418400 . 0 

0.139 

418400.0 ; 

CAE 

CAD 

4 


5 1 

0.139 

418400 . 0 

0.139 

418400.0 ; 

CAD 

CAC 

5 


6 1 

0.139 

418400 . 0 

0.139 

418400.0 ; 

CAC 

CAB 

6 


7 1 

0.139 

418400 . 0 

0.139 

418400.0 ; 

CAB 

CAA 

7 


8 1 

0.139 

418400 . 0 

0.139 

418400.0 ; 

CAA 

CAG 

8 


9 1 

0.133 

376560.0 

0.133 

376560.0 ; 

CAG 

NAH 

9 


10 1 

0.146 

251040 . 0 

0.146 

251040.0 ; 

NAH 

OAI 

10 


11 1 

0.100 

313800.0 

0.100 

313800.0 ; 

OAI 

HAA 


pairs 

] 





ai aj 

fu 

c0 

cl 



1 4 

1 


; 

CAJ 

CAD 

1 6 

1 


; 

CAJ 

CAB 

1 8 

1 


; 

CAJ 

CAG 

2 5 

1 


; 

CAF 

CAC 

2 9 

1 


; 

CAF 

NAH 


112 



3 

6 

1 




; CAE 

CAB 






3 

8 

1 




; CAE 

CAG 






4 

7 

1 




; CAD 

CAA 






5 

8 

1 




; CAC 

CAG 






6 

9 

1 




; CAB 

NAH 






7 

10 

1 




; CAA 

OAI 






8 

11 

1 




; CAG 

HAA 






angles 

] 












ai 

aj 

ak 

fu 

c0 


cl 








1 

2 

3 

1 

120.0 


418.4 120 

.0 

418.4 

; 

CAJ 

CAF 

CAE 


1 

2 

7 

1 

120.0 


418.4 120 

.0 

418.4 

; 

CAJ 

CAF 

CAA 


3 

2 

7 

1 

120.0 


418.4 120 

.0 

418.4 

; 

CAE 

CAF 

CAA 


2 

3 

4 

1 

120.0 


418.4 120 

.0 

418.4 

; 

CAF 

CAE 

CAD 


3 

4 

5 

1 

120.0 


418.4 120 

.0 

418.4 

; 

CAE 

CAD 

CAC 


4 

5 

6 

1 

120.0 


418.4 120 

.0 

418.4 

; 

CAD 

CAC 

CAB 


5 

6 

7 

1 

120.0 


418.4 120 

.0 

418.4 

; 

CAC 

CAB 

CAA 


2 

7 

6 

1 

120.0 


418.4 120 

.0 

418.4 

; 

CAF 

CAA 

CAB 


2 

7 

8 

1 

120.0 


418.4 120 

.0 

418.4 

; 

CAF 

CAA 

CAG 


6 

7 

8 

1 

120.0 


418.4 120 

.0 

418.4 

; 

CAB 

CAA 

CAG 


7 

8 

9 

1 

120.0 


418.4 120 

.0 

418.4 

; 

CAA 

CAG 

NAH 


8 

9 

10 

1 

109.5 


376.6 109 

.5 

376.6 

; 

CAG 

NAH 

OAI 


9 

10 

11 

1 

109.5 


397.5 109 

.5 

397.5 

; 

NAH 

OAI 

HAA 


dihedrals 

] 











ai 

aj 

ak 

al 

fu 

c0 

cl m c0 

cl m 






2 

7 

3 

1 

2 

0.0 

1673.6 0 

O 

o 

1673.6 0 ; 

IDI 

CAF 

CAA 

CAE 

CAJ 

7 

2 

8 

6 

2 

0.0 

1673.6 0 

o 

o 

1673.6 0 ; 

IDI 

CAA 

CAF 

CAG 

CAB 

2 

3 

4 

5 

2 

0.0 

1673.6 0 

o 

o 

1673.6 0 ; 

IDI 

CAF 

CAE 

CAD 

CAC 

3 

4 

5 

6 

2 

0.0 

1673.6 0 

o 

o 

1673.6 0 ; 

IDI 

CAE 

CAD 

CAC 

CAB 

4 

5 

6 

7 

2 

0.0 

1673.6 0 

o 

o 

1673.6 0 ; 

IDI 

CAD 

CAC 

CAB 

CAA 

5 

6 

7 

2 

2 

0.0 

1673.6 0 

o 

o 

1673.6 0 ; 

IDI 

CAC 

CAB 

CAA 

CAF 

6 

7 

2 

3 

2 

0.0 

1673.6 0 

o 

o 

1673.6 0 ; 

IDI 

CAB 

CAA 

CAF 

CAE 

7 

2 

3 

4 

2 

0.0 

1673.6 0 

o 

o 

1673.6 0 ; 

IDI 

CAA 

CAF 

CAE 

CAD 

2 

7 

8 

9 

1 180.0 

41.8 2 180 . 0 

41.8 2 ; 

Dl 

CAF 

CAA 

CAG 

NAH 

7 

8 

9 

10 

1 

0.0 

3.8 3 

o 

o 

3.8 3 ; 

Dl 

CAA 

CAG 

NAH 

OAI 

8 

9 

10 

11 

1 

0.0 

1.3 3 

o 

o 

1.3 3 ; 

Dl 

CAG 

NAH 

OAI 

HAA 


DZPanc.itp '. 

[ moleculetype ] 
; name nrexcl 
DZPANC 3 


This file was generated by PRODRG version 040913.0519 
PRODRG written/copyrighted by Daan van Aalten 

Questions/comments to dava@davapcl.bioch.dundee.ac.uk 

When using this software in a publication, cite: 

A. W. Schuettelkopf and D. M. F. van Aalten (2004) . 
PRODRG - a tool for high-throughput crystallography 
of protein-ligand complexes. 

Acta Crystallogr. D60, 1355--1363. 


[ atoms ] 


nr 

type 

resnr 

resid 

atom cgnr 

charge 

1 

CH3 

1 

DZPANC 

CAJ 

1 

-0 .017304 

2 

CB 

1 

DZPANC 

CAF 

1 

-0.0729 

3 

CR61 

1 

DZPANC 

CAE 

1 

-0.042 

4 

CR61 

1 

DZPANC 

CAD 

1 

-0.196 

5 

CR61 

1 

DZPANC 

CAC 

2 

0.0312 

6 

CR61 

1 

DZPANC 

CAB 

2 

-0.116 

7 

CB 

1 

DZPANC 

CAA 

2 

0 . 1421 

8 

CR61 

1 

DZPANC 

CAG 

2 

-0 . 1457 
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9 

NL 


1 

DZPANC 

NAH 

2 0.0353 






10 

OA 


1 

DZPANC 

OAI 

3 -0.6187 






bonds 












ai 

aj 

fu 

c0 


Cl 







1 

2 

1 

0.153 

334720.0 0.153 334720.0 ; 

CAJ 

CAF 




2 

3 

1 

0.139 

418400.0 0.139 418400.0 ; 

CAF 

CAE 




2 

7 

1 

0.139 

418400.0 0.139 418400.0 ; 

CAF 

CAA 




3 

4 

1 

0.139 

418400.0 0.139 418400.0 ; 

CAE 

CAD 




4 

5 

1 

0.139 

418400.0 0.139 418400.0 ; 

CAD 

CAC 




5 

6 

1 

0.139 

418400.0 0.139 418400.0 ; 

CAC 

CAB 




6 

7 

1 

0.139 

418400.0 0.139 418400.0 ; 

CAB 

CAA 




7 

8 

1 

0.139 

418400.0 0.139 418400.0 ; 

CAA 

CAG 




8 

9 

1 

0.133 

376560.0 0.133 376560.0 ; 

CAG 

NAH 




9 

10 

1 

0.146 

251040.0 0.146 251040.0 ; 

NAH 

OAI 




pairs 












ai 

aj 

fu 

c0 


cl 







1 

4 

1 




; 

CAJ CAD 






1 

6 

1 




; 

CAJ CAB 






1 

8 

1 




; 

CAJ CAG 






2 

5 

1 




; 

CAF CAC 






2 

9 

1 




; 

CAF NAH 






3 

6 

1 




; 

CAE CAB 






3 

8 

1 




; 

CAE CAG 






4 

7 

1 




; 

CAD CAA 






5 

8 

1 




; 

CAC CAG 






6 

9 

1 




; 

CAB NAH 






7 

10 

1 




; 

CAA OAI 






angles 

] 











ai 

aj 

ak 

fu 


c0 

Cl 







1 

2 

3 

1 

120.0 

418.4 

120.0 418.4 

; 

CAJ 

CAF 

CAE 


1 

2 

7 

1 

120.0 

418.4 

120.0 418.4 

; 

CAJ 

CAF 

CAA 


3 

2 

7 

1 

120.0 

418.4 

120.0 418.4 

; 

CAE 

CAF 

CAA 


2 

3 

4 

1 

120.0 

418.4 

120.0 418.4 

; 

CAF 

CAE 

CAD 


3 

4 

5 

1 

120.0 

418.4 

120.0 418.4 

; 

CAE 

CAD 

CAC 


4 

5 

6 

1 

120.0 

418.4 

120.0 418.4 

; 

CAD 

CAC 

CAB 


5 

6 

7 

1 

120.0 

418.4 

120.0 418.4 

; 

CAC 

CAB 

CAA 


2 

7 

6 

1 

120.0 

418.4 

120.0 418.4 

; 

CAF 

CAA 

CAB 


2 

7 

8 

1 

120.0 

418.4 

120.0 418.4 

; 

CAF 

CAA 

CAG 


6 

7 

8 

1 

120.0 

418.4 

120.0 418.4 

; 

CAB 

CAA 

CAG 


7 

8 

9 

1 

120.0 

418.4 

120.0 418.4 

; 

CAA 

CAG 

NAH 


8 

9 

10 

1 

109.5 

376.6 

109.5 376.6 

; 

CAG 

NAH 

OAI 


dihedrals 

] 










ai 

aj 

ak 

al 

fu c0 

cl m 

c0 cl m 






2 

7 

3 

1 


2 0.0 

1673.6 

0 0.0 1673.6 0 ; 

IDI 

CAF 

CAA 

CAE 

CAJ 

7 

2 

8 

6 


2 0.0 

1673.6 

0 0.0 1673.6 0 ; 

IDI 

CAA 

CAF 

CAG 

CAB 

2 

3 

4 

5 


2 0.0 

1673.6 

0 0.0 1673.6 0 ; 

IDI 

CAF 

CAE 

CAD 

CAC 

3 

4 

5 

6 


2 0.0 

1673.6 

0 0.0 1673.6 0 ; 

IDI 

CAE 

CAD 

CAC 

CAB 

4 

5 

6 

7 


2 0.0 

1673.6 

0 0.0 1673.6 0 ; 

IDI 

CAD 

CAC 

CAB 

CAA 

5 

6 

7 

2 


2 0.0 

1673.6 

0 0.0 1673.6 0 ; 

IDI 

CAC 

CAB 

CAA 

CAF 

6 

7 

2 

3 


2 0.0 

1673.6 

0 0.0 1673.6 0 ; 

IDI 

CAB 

CAA 

CAF 

CAE 

7 

2 

3 

4 


2 0.0 

1673.6 

0 0.0 1673.6 0 ; 

IDI 

CAA 

CAF 

CAE 

CAD 

2 

7 

8 

9 


1 180.0 

41.8 

2 180.0 41.8 2 ; 

Dl 

CAF 

CAA 

CAG 

NAH 

7 

8 

9 

10 


1 0.0 

3.8 

3 0.0 3.8 3 ; 

Dl 

CAA 

CAG 

NAH 

OAI 


6.2 APÊNDICE B: ADIÇÃO DA TOPOLOGIA DO SGA NO CAMPO DE FORÇAS 
GROMOS 87 


Para o 
neurotóxico 


GROMACS 3.2.1 reconhecer a ligação covalente entre o agente 
GA e o aminoácido Ser203 foi necessário a adição da topologia e das 
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coordenadas de um novo resíduo que chamamos de SGA (Ser 203 ligada ao GA) e, 
isto foi feito através do seguinte procedimento: 

1- Obtêm-se os arquivos com extensão itp e gro, pelo mesmo processo descrito 
no APÊNDICE A, para o resíduo SGA. 

2 - No arquivo aminoacids.dat, onde estão os nomes de todos os resíduos 
reconhecidos pelo campo de força GROMOS 87, incluído no diretório onde o 
GROMACS 3.2.1 foi instalado, adiciona-se o nome para um novo resíduo, 
SGA, no caso deste trabalho. 

3 - Copiam-se as topologias referentes a nome de átomos, cargas atômicas já 
refinadas, ligações, ângulos, torcionais próprios e torcionais impróprios, 
presentes no arquivo de extensão itp, para dentro do arquivo ffgmx.rtp, 
seguindo o modelo de outros já incluídos, onde estão todos os parâmetros 
reconhecidos pelo campo de forças GROMOS 87. 

4- Disto tudo, o resíduo SGA passará a ser reconhecido pelo GROMACS 3.2.1 . 

Parte do conteúdo do arquivo ffgmx.rtp, onde foram adicionadas as topologias 

do SGA, é mostrado a seguir: 

[ SGA ] 

; ELIMINEI O HAA E O HAC E CHAMEI O HAB DE H 
[ atoms ] 


C6 

CH3 

-0 

. 0805 

0 

C5 

CH2 

0 

. 4248 

0 

04 

OS 

-0 

.5295 

0 

P 

P 

0 

.7735 

1 

NI 

NL 

-0 

.7584 

1 

C4 

CH3 

0 

. 1872 

2 

C3 

CH3 

0 

.8222 

2 

03 

OM 

-0 

.8328 

2 

OG 

OS 

-1 

. 0859 

2 

CB 

CH2 

1 

. 7243 

2 

CA 

CHI 

-0 

.4609 

3 

C 

C 

1 

. 1769 

3 

0 

0 

-0 

. 9633 

3 

N 

N 

-0 

.5946 

3 

H 

H 

0 

.1970 

3 


bonds 

] 




C6 

C5 

; 0 . 

.153 334720.0 0. 

.153 3347 

C5 

04 

; 0 . 

.143 251040.0 0. 

.143 2510 

04 

P 

; 0.1 61 

251040.0 0.161 

251040 . 0 

P 

NI 

; 0.1 62 

304571.8 0.162 

304571 . 8 

P 

03 

; 0 . 148 

376560.0 0.148 

376560.0 

P 

OG 

; 0.161 

251040.0 0.161 

251040 . 0 

NI 

C4 

; 0 . 147 

376560.0 0.147 

376560.0 

NI 

C3 

; 0 . 147 

376560.0 0.147 

376560.0 

OG 

CB 

; 0 . 143 

251040.0 0.143 

251040 . 0 

CB 

CA 

; 0 . 153 

334720.0 0.153 

334720 . 0 

CA 

C 

; 0 . 153 

334720.0 0.153 

334720 . 0 

CA 

N 

; 0 . 147 

376560.0 0.147 

376560.0 
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C 0 ; 0 . 123 502080.0 0.123 502080.0 

N H ; 0 . 1 0 0 374468.0 0.100 374468.0 


-C N 
[ angles ] 


C6 

C5 


04 


; 111 . 0 


460.2 

111 . 0 


460.2 

C5 

04 


P 


; 120 . 0 


397.5 

120 . 0 


397.5 

04 

P 


NI 


; 103.0 


397.5 

103.0 


397.5 

04 

P 


03 


; 1 0 9 . 6 


397.5 

109.6 


397.5 

NI 

P 


03 


; 1 0 9 . 6 


397.5 

109.6 


397.5 

04 

P 


OG 


; 103.0 


397.5 

103.0 


397.5 

NI 

P 


OG 


; 103.0 


397.5 

103.0 


397.5 

03 

P 


OG 


; 1 0 9 . 6 


397.5 

109.6 


397.5 

P 

NI 


C4 


; 1 0 9 . 5 


376.6 

109.5 


376.6 

P 

NI 


C3 


; 1 0 9 . 5 


376.6 

109.5 


376.6 

C4 

NI 


C3 


; 1 0 9 . 5 


376.6 

109.5 


376.6 

P 

OG 


CB 


; 120 . 0 


397.5 

120 . 0 


397.5 

OG 

CB 


CA 


; 111 . 0 


460.2 

111 . 0 


460.2 

CB 

CA 


C 


; 1 0 9 . 5 


460.2 

109.5 


460.2 

CB 

CA 


N 


; 1 0 9 . 5 


460.2 

109.5 


460.2 

C 

CA 


N 


; 1 0 9 . 5 


460.2 

109.5 


460.2 

CA 

C 


O 


; 121 . 0 


502.1 

121.0 


502 . 1 

CA 

N 


H 


; 1 0 9 . 5 


376.6 

109.5 


376.6 

[ dihedrals 

] 








P 

04 

C5 

C6 

; 0 . 0 

3.8 

3 0 

.0 3.8 

3 


C5 

04 

P 


OG 

; 0 . 0 

1 . 0 

3 0 

.0 1.0 

3 


C5 

04 

P 


OG 

; 0 . 0 

3.1 

2 0 

.0 3.1 

2 


04 

P 

NI 

C3 

; 0 . 0 

1 . 0 

3 0 

O 

h- 1 

o 

3 


04 

P 

NI 

C3 

; 0 . 0 

3.1 

2 0 

.0 3.1 

2 


CB 

OG 

P 


04 

; 0 . 0 

1 . 0 

3 0 

o 

h- 1 

o 

3 


CB 

OG 

P 


04 

; 0 . 0 

3.1 

2 0 

.0 3.1 

2 


P 

OG 

CB 

CA 

; 0 . 0 

3.8 

3 0 

.0 3.8 

3 


OG 

CB 

CA 

N 

; 0 . 0 

5.9 

3 0 

.0 5.9 

3 


CB 

CA 

C 


O 

; 0 . 0 

0 . 4 

6 0 

.0 0.4 

6 


[ impropers 

] 








P 

04 

03 

NI 

; 3 5 . 3 

836.8 

0 35 

.3 836.8 

0 


NI 

P 

C3 

C4 

; 3 5 . 3 

836.8 

0 35 

.3 836.8 

0 


CA 

CB 

N 


C 

; 3 5 . 3 

836.8 

0 35 

.3 836.8 

0 



6.3 APÊNDICE C: PROGRAMAS DO PACOTE GROMACS 3.2.1 

A seguir são apresentados os programas executáveis do pacote GROMACS 
3.2.1 utilizados para gerar os arquivos de entrada, executar os cálculos de 
minimização de energia e de dinâmica molecular e extrair informações dos arquivos 
de saída para posterior análise. 

C.1 CONSTRUÇÃO DOS ARQUIVOS ENTRADA 
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Para a construção dos arquivos entrada, foi necessário somente o arquivo 
HuAChE+GA.pdb, contendo as coordenadas tridimensionais da proteína com o 
resíduo GA ligado covalentemente ao aminoácido Ser203. A seqüência de 
programas utilizados foi a seguinte: 

pdb2gmx -f HuAChE+SGA . pdb -o sistema. gro -p sistema. top 
editconf -bt cubic -f sistema. gro -o sistema. gro -c -d 0,65 
genbox -cp sistema. gro -cs spc216.gro -o sistema_st . gro -p sistema. top 
make_ndx -f sistema. tpr 

O programa “pdb2gmx” converte o arquivo .pdb da proteína no arquivo 
sistema. gro, contendo as coordenadas da proteína e do SGA no formato que o 
pacote GROMACS 3.2.1 reconhece, e no arquivo sistema.top contendo os dados de 
topologia da HuAChE+GA ; 

O programa “editconf” gera em torno da proteína, já com o resíduo SGA incluído, 
uma caixa de tamanho adequado ao volume da mesma; 

O programa “genbox” solvata o sistema, isto é, introduz nesta caixa as moléculas 
do solvente; 

O programa “make_ndx” é utilizado para criar os grupos dentro do sistema para 
posterior análise ex: grupo proteína (já com o SGA incluído)/2PA/W. 


C.1.1 Incluindo os Arquivos Correspondentes ao 2pam, Dzp e Dzpanc no 
Sistema Huache+Sga 

Os arquivos 2pam.itp, 2pam.gro, dzp.itp, dzp. gro, dzpanc. itp e dzpanc. gro 
gerados pelo servidor PRODRG e já com as cargas corrigidas por métodos 
quânticos, foram incluídos no sistema segundo o procedimento: 

1- Editar os arquivos sistema. gro, 2pam/dzp/dzpanc.gro. 

2- Copiar e colar as coordenadas da 2pam, dzp e dzpanc no final do arquivo 
sistema. gro. 
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3- Executar editconf para acertar a numeração. 

4- Corrigir o número total de átomos no topo do arquivo sistema. gro. 

5- Editar o arquivo sistema.top. 

6- Querendo-se simular, por exemplo, o sistema com a 2PAM, acrescentar 
“#include 2pam.itp” abaixo da linha “Include topologies”. 

7- Acrescentar na última linha do arquivo os nomes e a quantidade de 
moléculas adicionadas de cada espécie (o mesmo nome deve constar do 
arquivo de topologia). 

C.1.2 Inserindo íons para Equilibrar a Carga do Sistema 

Como os sistemas estudados não apresentaram carga líquida zero, foi 
necessário incluir íons para equilibrar as cargas. Para a HuAChE, que apresentou 
cargas negativas em todos os casos, foram adicionados íons Na + . Estes íons foram 
inseridos no sistema segundo o procedimento descrito abaixo: 

1- Editar o arquivo de topologias, sistema.top. 

2- Digitar o nome e o número de íons que se deseja colocar, abaixo da linha 
“[Molecules]”. 

3- Executar o programa genion para gerar os íons. 

genion -s sistema_st . tpr -o sistema_st . gro -pot -nice 0 -nn "número de ions" 
-nname "simbolo quimico do ion" -nq "carga do ion" 

4- Na opção “Select a group”, selecionar sempre o grupo correspondente ao 
solvente. 

5- Subtrair no arquivo .top, abaixo da linha “[Molecules]”, o número de 
moléculas do ion adicionado do número de moléculas do solvente. 

C.2 MINIMIZAÇÃO DE ENERGIA E DINÂMICA MOLECULAR 

Uma vez construídos os arquivos de entrada, a próxima etapa será a execução 
dos programas para a minimização de energia e simulação de dinâmica molecular. 
A seqüência de programas utilizados foi a seguinte: 

Para a minimização de energia : 

grompp -s st.mdp -c sistema_st . gro -p sistema.top -o sistema_st . tpr 

mdrun -v -s sistema_st . tpr -o sistema_st . trr -c sistema_st . gro -g stlog>&st . j ob& 

Para dinâmica das moléculas de água : 

grompp -v -f pr.mdp -c sistema_st . gro -p sistema.top -o sistema_pr . tpr 
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mdrun -v -s sistema_pr . tpr -o sistema_pr . trr -c sistema_pr . gro -g prlog>&pr . j og& 

Para a dinâmica molecular : 

grompp -v -f md.mdp -c sistema_pr . gro -p sistema. top -o sistema_md . tpr 

mdrun -v -s sistema_md . tpr -o sistema_md . trr -c sistema_md . gro -g mdlog>&md . j ob& 

O programa grompp concatena dados dos arquivos de coordenadas (.gro), de 
topologia (.top) e de parâmetros (.mdp) do sistema em um só arquivo que servirá de 
arquivo de entrada para a execução dos cálculos que se seguirão. 

O programa mdrun dá início as simulações e gera os arquivos de saída .trr e . 
edr, que contém todos os resultados dos cálculos, os arquivos .gro com as 
coordenadas dos últimos quadros gerados, nas respectivas etapas e os arquivo de 
st/pr/mdlog e Job que contém dados referentes ao andamento dos cálculos. 

São os arquivos .mdp que contém os parâmetros do cálculo a ser realizado e, 
para cada etapa do cálculo, foi utilizado um arquivo .mdp diferente. No item C.4 são 
apresentados os conteúdos dos arquivos .mdp utilizados neste trabalho. 

C.3 PROCEDIMENTO PARA A ANÁLISE DE RESULTADOS DAS DINÂMICAS 

Depois de concluída a dinâmica molecular, foram executados os programas para 
gerar os gráficos de variação de energia, DRMQ temporal e espacial, número de 
ligações hidrogênio, distâncias interatômicas e para extrair quadros ao longo da 
dinâmica. Estes programas são apresentados abaixo: 

Para gerar o gráfico de variação de energia total : 

g_energy -f ener.edr -s sistema_md . tpr -o energia. xvg 

escolhe-se a opção 12, opção referente ao cálculo de energia total. 

Para extrair quadros dos sistemas HuAChE- GA com a 2PAM, DZP ou a DZPanc : 

make_ndx -f sistema_md . tpr [usar o h (help) para criar o grupo 2PAM, DZP ou 
DZPanc com o sistema, por exemplo: 2PAM/sistema ]. 

trjconv -fit Progressive -n index.ndx -f sistema md.trr (ou traj.xtc) -s 
sistema_md. tpr -o sistema. pdb -b t 0 -e t f 

Em que t 0 é o quadro no tempo t 0 ps e t f o quadro no tempo t f ps, ou seja, 
extraem-se as coordenadas indo do quadro t 0 a t f ps. 

A opção -fit Progressive ajusta os quadros progressivamente ao longo da 
dinâmica. Isso evita distorções que podem ser provocadas pelo possível 
deslocamento do sistema dentro da caixa d’água. 
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Para gerar os gráficos de DRMQ temporal e DRMQ espacial : 

Cria-se primeiro os grupos que se quer analisar com make_ndx\ 

make_ndx -f sistema_md . tpr 

DRMQ temporal : 

g_rms -n index.ndx -s sistema_md . tpr -f sistema_md . trr -o DRMQtemp.xvg -prev 1 

DRMQ espacial : 

g_rmsf -n index.ndx -s sistema_md . tpr -f sistema_md . trr -o DRMQesp.xvg -res 

Para gerar os gráficos de distâncias interatômicas : 

Criar grupos de pares de átomos com make_ndx: 

make_ndx -f sistema_md . tpr 

Selecionam-se os átomos da seguinte forma: 
an <ENTER>, 

para tantos átomos quantos forem necessários para a análise, onde n é o 
número do átomo que pode ser encontrado no arquivo sistema_md.gro. 

Executar o programa g_dist 

g_dist -f traj.xtc -s sistema_md . tpr -n index.ndx -o distancia . xvg 

Escolhe-se então, os pares criados anteriormente pelo make_ndx para cada 
distância interatômica a ser calculada. 

Analisa-se o arquivo .xvg com o programa XMGRACE. 

Para gerar os gráficos de número de ligações hidrogênio : 

Cria-se os grupos a serem analisados com make_ndx\ 

make_ndx -f sistema_md . tpr 

Executa-se o programa g_hbond da seguinte forma: 

g_hbond -f traj.xtc -s sistema_md . tpr -n index.ndx -num NumLigH.xvg , onde 

NumLigH.xvg é o arquivo que contem o número de ligações hidrogênio formadas em 
função do tempo de simulação. 

Todos os arquivos com extensão .xvg contidos neste trabalho, foram visualmente 
analisados pelo programa XMGRACE e melhorados no programa CALC, incluso no 
pacote OPENOFFICE. 

C.4 CONTEÚDO DOS ARQUIVOS DE PARÂMETROS .MDP 

st.mdp : 

; User ***** (236) 

; Wed Nov 13 17:12:44 2004 
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; Input file 

title 

cpp 

define 

constraints 

integrator 

dt 

nsteps 

nstlist 

ns_type 

rlist 

coulombtype 

rcoulomb 

rvdw 

f ourierspacing 
fourier nx 
fourier ny 
fourier nz 
pme_order 
ewald rtol 
optimize fft 


PEP 

/lib/cpp 

-DFLEXIBLE 

none 

steep 

0.01 ; ps ! 

20000 
10 

grid 
0 . 9 
PME 
0 . 9 
0 . 9 
0 . 12 
0 
0 
0 
4 

le-5 

yes 


; Energy minimizing stuff 

r 

emtol = 418.4 

emstep = 0.01 


pr.mdp : 


ps ! 

total 50 ps. 


; User ***** (236) 

; Wed Nov 13 17:12:44 2004 

; Input file 

title = Dinamica 

cpp = /lib/cpp 

define = -DPOSRES 

constraints = all-bonds 

integrator = md 

dt = 0.002 

nsteps = 25000 ; 

nstcomm = 1 

nstxout = 2500 

nstvout = 5000 

nstfout = 0 

nstlog = 5000 

nstenergy = 5000 

nstlist = 10 

ns_type = grid 

rlist =1.5 

rcoulomb = 1.5 

rvdw = 1.5 

; Berendsen temperature coupling is on in two groups 
Tcoupl = berendsen 


tc-grps = Protein * 

tau_t = 0.1 

ref_t = 300 

; Energy monitoring 
energygrps = Protein * 

; Pressure coupling is not on 

Pcoupl = no 

tau_p = 0.5 

compressibility = 4.5e-5 

ref_p = 1.0 


SOL 


Na 


0 . 1 


0 . 1 


0 . 1 


300 


SOL 


300 


Na 


300 
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is on at 300 K. 
yes 
300 . 0 
173529 


md.mdp : 


; Generate velocites 
gen vel 
gen temp 
gen^seed 


User ***** (236) 

Wed Nov 13 17:12:44 2004 
Input file 


title 

cpp 

define 

constraints 

integrator 

tinit 

dt 

nsteps 

nstcomm 

nstxout 

nstvout 

nstf out 

nstlog 

nstenergy 

nstxtcout 

nstlist 

ns_type 

rlist 

coulombtype 

rcoulomb 

rvdw 

f ourierspacing 


ache MD 
/lib/cpp 
-DFLEX_SPC 
none 
md 
0 

0 . 002 
500000 
1 

10000 
10000 
0 

100 
100 
500 
5 

grid 
1 . 5 
PME 
1 . 5 
5 


ps ! 

total lOOOps 


= 1 

= 0.12 


fourier nx 

= 0 




fourier ny 

= 0 




fourier nz 

= 0 




pme order 

= 4 




ewald rtol 

= le-5 




optimize fft 

= yes 




; Berendsen 

temperature coupling 

is on 

in three 

groups 

Tcoupl 

= Berendsen 




tau t 

= 0.1 

0.1 

0 . 1 

0 . 1 

tc-grps 

= protein 

★ ★ ★ 

SOL 

Na 

ref t 

= 300 

300 

300 

300 


Pressure coupling is not on 


Pcoupl 
pcoupltype 
tau_p 

compres sibility 
ref_p 

; Generate velocities is on at 300 K 
gen vel 
gen temp 
gen_seed 


= no 

= isotropic 
= 2 

= 4 . 5e-5 

= 1.0 


yes 
300 . 0 
173529 


Em que os asteriscos vermelhos podem ser a 2PAM, DZP ou a DZPanc, 
dependendo da simulação. 
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